Mathematies 
and 
Mathematical 
Physics 








SECTION B 


JOURNAL OF 
RESEARCH 


NATIONAL BUREAU OF STANDARDS 








Journal of Research 
of the 
National Bureau of Standards 





Editor: Chester H. Page 


Associate Editors: F. L. Alt, E. L. Crow 
Publication dates: Stan, 35 ne 10, Sept. 10, Dec. 10, 





. 


| 


papers of interest i oo | to 
times pred bscripti bread 
a % su ion: c, 
$4.00; foreign, $4.75. 


mainly fet the mathematician ond the theo 


subscription: domestic, $2 S tore, $2.73 


Functions and Activities 

The functions of the National Bureau of Standards are set forth in the Act of 
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Physical systems and experimental situations ars entitv g must have a significan t is independent 
quantit itively described in terms ol observables ol our choice of Measurement unit The observable 
entities Which are s ibyject to mMensure ent I} Ost must be ivarinnt to our mode ol aescription If we 
entities are of various kinds, e.g., force, mass, charg change our Measurement unit from g,. to ¢ where 
ete wo entities are said to be of the same kind 
if thev can be compared Lt if thre q lestiol / / s a times as ls ‘ 
is Meanimetul For example, we can compare oné 
foree wit! nother to determine whicl is large! Dut le ; —— Po a — 
we cnnnot sav whethel ivivenh distance Ss iurger ol 
smaller than i particular mass _ ; 
The mplicatior Ol thre above remarks Ss that qj qj r 
phvsicat entitles Cal be sorted into ordered ¢| Ws 
or one-dimensional sets, of like entities where q denotes the measure of the same sample 
of g relative to g,. as Measurement unit 
¢ y Bi experimental laws of physics are expressibl 
s proportions among Measures OF entitles such as 
i 
he entities i ire snid to be of t ve od! 4 
| il cs ciilterel ( ~ not be com 
pared of different | ais Che commor ‘ proportionality constant Cle pe ds ipon the 
propel of « tities in the g-chai , CO nic thy CHOICE ot menusurement ubtt iO ey ous entities 
he called ‘‘q-ness’’. for an identifving label nvolved ind upon the configuration of the experi 
Obs byle 1 ot onl ibie » ord le 1 st p It is customal >. 7 nk of the 
— bu sans lso subiect to measurement rhus eusures on t right as rel to more “‘elemen 
vs EXD mentally detern p 1] o hetwes entities Whose measu nent nits ive 
vo ent , \ ‘ Uready been chosen The measurement unit of the 
. new titv ¢g is thel lenne ) Ssicl con 
e ‘ , re ven t Val to f ol is rhe ead eCXp ( ta 
| itiol 
The experimental procedures of physical science are \s we s HT see later, a few) surement units can 
Dnased Oo | ( rpbitral ehoie ota particulal sample bn writ I L\ ehosen is L Dasis [ol 1 sel ol nits 
of an entitv as a easurement u for entities of that These basic units, together with assigned values for 
land Why eneral sample of that na s the the &’s in the fundamental equations, determine the 
described bv the state { remaining measurement nits ermanhe to the 
svstem 
er oe eee It should be noted that, during the development 
ol il SCLTCTICS thre factol } nel (Les thie effects ol 
- symbol fo e partic rsample chosen | ¢Xperimental variables that have eit er inadvertently 
is called the value, or measure, of the sample, g the interaction forces among current-cal ving con 
In contrast. the Moh hardness seale is a familiar | @uctors depend upon the fluid medium filling the 
example of chain ordering vthout metriu properties laboratory he permeability ot an varies so little 
lf ne concept of al observabl , has any real with normal environment il changes, that the effect 
’ or is to have any quantitative usefulness, the of this variable would not appear in exp riments 
~ leading to an initial formulation of Ampere’s force 
law 








It should also be noted that there are hidden 
dependencies among the k’s. For example, the in- 
teraction force between two stationary charges, the 
interaction force between two current loops, and the 
relation between current and charge implied by the 
magnetic effect of a moving charge, are not inde- 
pendent. The three k’s involved satisfy a relation 
involving the speed of light in the medium in which 
the experiments are performed. The existence of a 
relation between a set of k’s and the measure of some 
physical entity implies that at least one & contains a 
hidden reference to a physical entity not explicitly 
referred to in the experimental measure equations 

There is a uniqueness postulate in physics, usually 
not stated, that 1s relevant to the completeness of a 
law. “If ig} depends upon ‘p), «s', ete., then any 
entity whose measure depends upon the same varia- 
bles in the same way, and upon no others, is an 
entity of the same kind as g."" Any apparent viola- 
tion of this postulate is an indication of hidden 
“variables” in the proportionality constants. These 
“variables” may be properties of nature which are 


beyond our power to change, such as the speed of 


light in vacuo. This constancy prevents the experi- 
menter from including these variables in any particu- 
lar law, and also allows the theoretician to distribute 
the effect of hidden. variables in any convenient 
self-consistent manner. 

The development of physics as a science, rather 
than as a taxonomy of empirical relations, was due 
to the abstract representation of these relations by 
mathematical relations among mathematical ele- 
ments, subject to manipulation by algebraic and 
differential processes. Maxwell wrote: ? 

“There are two methods of interpreting the equa- 
tions relating to geometry and other concrete sciences. 

“We may regard the symbols which occur in the 
equation as of themselves denoting lines, masses 
times, &e.; or we may consider each symbol as 
denoting only the numerical value of the correspond- 
ing quantity, the concrete unit to which it is referred 
. being tacitly understood 

“If we adopt the first method we shall often have 
difficulty in interpreting terms which make their 
appearance during the calculations. We shall there- 
fore consider all the written symbols as mere numeri- 
eal quantities, and therefore subject to all the opera- 
tions of arithmetic during the process of calculation 
But in the original equations and the final equations, 
in which every term has to be interpreted in a physi- 
cal sense, we must convert every numerical expression 
into a concrete quantity by multiplying it by the 
unit of that kind of quantity.” 

We next examine the basis of the representation 
of physical entities (concrete quantities) by mathe- 
matical elements, and the structure of the associated 
mathematics. For simplicity, we shall consider all 
physical entities as scalars; 
be so described in terms of their components 

Since samples of q can be chain-ordered, these 
samples can be represented by points in a one- 


vectors and tensors can 


dimensional topological space. In addition, the 
? Encyclopaedia Britannica, Ninth Edition, VII, 24 
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measurability of observables implies that this is a 
metric space; and in particular, the measurability 
as a ratio to a Measurement unit implies au one 
dimensional abstract vector space We use the same 
svmbol q to denote both a physical entity and its 
representative point in the vector space. In its 
latter role as a mathematical element, the symbolic 
quantity g is often ealled a “physical quantity 

This is the reason I have been using the term ‘‘e1 
tity’ for the physical Cuse, 
an “entity” and its representation by a “quantity 

The duality of the svmbol q causes little difficulty 
the ambiguity is usually 
On the other hand, as we shall see later, substitution 


to distinguish between 


resolved by the context 


of a name for q can cause contusion, such as the Ist 
ol “ampere” for the name of the measurement unit 
of current as well as for a mathematical element 


The particular ent ty sam ple » Bas used for a measu 


ment unit, will he represented hy a point denoted by 
q), the hase or unit ctor of the Space i] wh ch thre 
entity is represented The vectol q Is al abstrae 


or symbolic unit, as contrasted with the more cor 
crete Measurement unit Soth abstract and meas 
urement units are often referred to as “units’’, with 
no modifier. Thus, to an experimenter, the unit o 
current is the ampere, and is specified by operattol ul 
procedures, while to Tans 
current is also the ampere, but is defined to be 
mathematical element! 

The relation between an entity 
measurement unit, “g Is |g 


theorists, the unit o 


sample and its 


tines as large us ¢ is 


/ 


represented by the equation 


¢ Gi\¢ 


This ts an equation relating points or veetors in the 


space of g-elements Let {[g] be a symbol for this 
q-space, the representation space for entities of thi 
kind g We can think of [g] as a symbol for the 
attribute of ‘“g-ness”’ of the entity 


Later q| Will 
also be called the ‘““dimension”’ of q 
Generalization of Measure Equations 


If k contains no hidden “variables’’, the theo 





retician generalizes the experimenter’s measur 
equation, 7 } Pp s+, to the abstract relatior 
q=kps, where q, p, and s are now mathematical 
elements in the lq], (pl, and [xs spaces The innocent 
looking “product’”’, ps, is not vet defined. Howeve1 


our postulate on the completeness of a physical law 
in combination with the assumed relation between 
qd. Ps and S, implies that the product operation has 
the following propert\ 

“The product of any point 
is a point in [g], and conversely, any point in 
can be expressed as the product of some point im [p 
and some point in [s]."". We can therefore defin« 

to be the product of the spaces [p] and [s] 


Nn p| and any point in 


ls 





q\=|p) Is}. 


Since the order ol changing the Pp and s conditions i 
an experument Is immaterial, we can require Ut! 


ibove 


 COMMULALIVE called ohe ent set bsence of coeflicients 
m rhe result Gg | I equ 0 shows that the 
ol nt wu Ss form a group, wit! | is thi identity 
lis u \ olsidered as the unit vector of the space 
Note that we are not here defining the result |g] of ol numerics This make the set ol ymbolic 
the product operation, but are defining the product units coherent With the set ol corresponding 
operatiol bv its result Space \ set of germal measurement units 1s 
We ilso define the less Ireque! tly needed sum represent d by a set of cohere t symbol nits 
of spaces to be ogical sum Since [p] is the spac Wh we treat physics from the symbolic, o1 
of p-ness, ic., the set of all possible values of p | quantity calculus, point of view, we consider the 
collect vel ( be defined to be the set of quantity eg lations as bast ind the measure equa- 
ill eleme! KOSSECSSI! -NeSs plus tiie set ol all Lions is derived The qual | equations unré 
eleme! S possessil PSS * Lf the set ol ll elemet ts theretore nval nt to choice of u ts "| Is 
DOSSCSS ( Cs ol ! ss ol bot! lt 
follow | / }: , 
( vstems 1, 2, a re different Kor 
vamp 
( mice Oo mb oO Ni ol bore 0 Thhe hewton 
ler ) e equ measure Of mass po cis) tine pound 
(Also time measure of acceleration mn centimeters per 
econa hes ct Linnie ‘ ded by second 
Sine element of multiplied b squi | 
‘ ! 0 \ ( i element I 
| sel onsistency oO ( experimental propo! Ss equation cd i be tactored 0 
onalities of plivsies unphes a unique multiplication 
by] hor titre ilo qdihbtact t! t tire iorm [ j ( 
JTOUDp \ I is Its 1entityv elem 
| 0 eral I may nvolve hidader ‘variables f , 
at mat il model must have ternal consist 
ev, and to be useful, should satisfy a e lence aS ind \a - ats 
. } = 
( ( Phi measu _ equation aeaduced irom i / / VBI n ‘ a 
mbo equattol to be identical with t! exper 
i measure equation being modeled The ind the conventional measure equatiol } 
sure equation kip s generalized to htel. seeks 
, / t / l 
( [kp 
7 Pe Se CY In fact, the invariance of the s bolic quantities 
ere $a unit vector ibitrary space. | i» and a yields directly 
Ce 4 ‘ we immediatel find the init 
al avn wna ! a 
() | ) 
so that / NH a san cone init equation, 
space relatio he above noncoherent units eld the measure 
equ on 
| . r m}2{a}a/a 
me 10gK —— for introducing various [1 In textbooks, the braces and subscripts are omitted 
es Was no WW dely recognized ul | tiie MIKSA | 
nd replaced by marginal notes 
em o mits Dbecume populal hae ad mensional 
VSis procedures ommonly used by \IKSA / 10 lp/idt- F volts. él maxwells 
T es brought to lel the hidden ‘variables’ 
e cony tio! il proportionality constants I} ! eoherent set ol mits, the equation has the 
const on k's, such as kA s general simple form 
to ke ip f 
KT) eo I 
il 


pond r « 


onstramt on thre issighinent of 
the spaces and the unit vectors of the T's ut 
troduced for completeness of the theory 
t of svmbol mits which makes thr Lheoretl il 
perimental measure equations identical is 


l‘o LTTLTLA TZ 


follow ing 


Lhe product 


show the 


| eps defines the prod wt ol p and 
2 ] p ‘ yields the multiplication table of 
tiv group of spaces O1 dimensions 
ol its equivalent 
i p vields the multiplication table 
for the group of coherent units 








Proportionality Constants implies a change of some geometric unit. For ex 
ample, the magnetizing force at the center of 


We have seen that the theoretician must use A(T), | eyrrent circle is 
or an ae op where the experimenter uses a H—kI) 
numeric The conventional way of including the 
hidden “variable’’ is to make ka nonnumeric. This | with the rationalized and nontationalized measures 
is done by “correcting” the experimenter’s equation 
to k},=1/2 
q} k > 18 
7 / L on 
and using the theoretical equations: 
to have fixed meanings, 


If H and k are 


q kps 
(q)=(kp)Xs). MH Wil 
. . ' , k} tk kh ' 
A given experimental law, such as Coulomb’s Law | viajdine : 
for electric charges, . . 5 
kh hk , 
FP} cig iq}/ir jh j 
is frequently subjected to different choices of pro- Coherence demands 
portionality constant by different writers. To H bd] 
correlate all these choices, and make them yield the :' 
same Meaning, requires interpreting them as the H AWAAYE 
measure equations relevant to a single symbolic | go that 
equation, with different choices of units. Thus H H I: 
bir 
F=ko.qo/r2 Hi H hk 
ields an 
yiele I rp? I q q Che relation k tor I suggests that kis a solid : 


and angle, measured in terms ol spheres or steradians 
us units. 
Perucea ‘emphasized that the factor 47 appearing 


in the ith system of coherent units. We assume a | ! nonrationalized equations of electromagnetism 


common choice of units for F and r. Since the left | “of &@ numeric, but the quantity 47 steradia 
hand sides of these equations are independent of /, Young ® introduced a proportionality constant, S 
the choice of unit system, the right hand sides must 
also be invariant with the choice. Therefore, 


representing the measure ot a sphere mto Various 
equations. 
Chis interpretation of rationalization is based o 


q he a fired symbolie equation, with the measures ex 
qr \ k ¥ pressed relative to rationalized o7 nonrational (| ‘ 
units If, on the other hand, we consider rational! 
which agrees with the direct interpretation, by an | zation to be a change of the symbolic equation 
experimenter, of the change of proportionality con- 
stant. But we also have k=({k},(k),, which is H,=2nI/ 
invariant, so that the change of the unit of g implies H.=1/2) 
a change of the (symbolic) unit of k. (Since (#) was | \., — - 
introduced for theoretical convenience, we do not : fi trl 
insist on the existence of a measurement unit, / i 
corresponding to the symbolic unit k -/ so that I] and HT clo hot represent the same ice ntu 
We can interpret the above as follows. The entity. If the factor 47 is interpreted as a geo 
germane Measurement unit of an entity (g) depends | metrical quantity, //7, and //, are not strictly of thi 
upon the choice of proportionality constant ({/}) in | same kind (dimension We can still satisfy the 
the defining equation for the germane unit. If the | perimenter, and= avoid paradoxes, by arbitra 
symbolic quantity g is to have the same meaning in | getting 
all cases, there must be a corresponding change in H H).. {H. I]. 
the unit of &, i1.e., in the unit of some real or fictitious 
entity upon which & depends. so that the experimenter need not distinguish bi 
It has earlier been mentioned that the experimental tween his rationalized measure of 7 and the measur 
proportionality constant, {k , depends upon the geo- | of a new “rationalized /7” The problem will be 
metrical configuration of the defining experiment. | discussed in detail in a later section. 
When & contains no hidden “variables” in the usual aan 
physical sense, we can interpret it to contain a me ghee 4+ se stent MESA t : ths 
geometrical quantity, represented by some space [/] E. Perucca, Solid angle in electromagnetist \ 
aa : . : quantities, La Ricerea Scientifiea 30, Suppl. 12, p. 220 il 
lhe change of (£) associated with a change of {k} now | “Leo Young, Elect IEE 75. ti 
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Units and Dimensions logarithms, angles, and circular and elliptic fune- 
; tions \laxwell, nowevel ilso said low cit 
\ given stage im oul knowledge of paysics, “By know the dimensions ol Vv quantity we are 
we have independent proportionalities involving able at once to deduce its numeri os wales a6 am 
\V kinds of entity The proportiol ilities iunply pressed terms of one system of units from its 
= vedi | numer lue as given in terms ol othersystem 
We ea Noose i is arbitrary numbers La . lac parap rased t They ¢ eat if 
dN irbitrary measurement units CV-n prototype for toll ao Actas the emankeal e of 6 quantita 
standards These nits are called “basi es ee eee <<  Beeke euall saammenenl iim 
they serve as basis of a unit system whose othe! sh naiad fo anessethed che And later (p. 6 
on © unite lt ved from these \ One of “The wi oft temperature may D ssigned mde 
thi probl sola sl indards laboratory top echelo pendently ol thre units of the entitres of mechanics 
s to build experimental set-ups for ; lization Consequently, the svmbol [7 denotu tempera- 
ort di ed units ture is rr ded as one ot the damental dimen- 
We ar pt lo THINK O ' tet log ee sions The dimension of ry riabl is 
cle Kelvi nd candela as the only ts defined | qanoted by 
Dy prototype § caras The mcasurenn it for Sie I measure of al e aepends upol Lhe 
ele s, however, det ed DS pro type stand ehore of untt e.g , degree o qhia d the chores 
Pa C CICK > evors define Uh cle to be of unit 1s dependent of the choices made for the 
i real on of 360 sof angle: ser sts define atin teasin eaalin. Gh eomeieale : ol = ial 
ade u on of 2m ° Simi Vy, sol I ess tially numer il put e dimensions |¢ 
: 7 spneres Bt a 2 = solid ad i822) o | ie i orl d solu ( 0 dinarily 
— ' suppressed o1 ored 
Not we made .\ ~atnel “bee : ll In appendix 2, 1t Is shown that [@ sin 6] 6 
h \ prototyp | — \ ; —" ne cos O I so that 4] lj is al illowabl ussivgh 
split ip ' tly For Cx upl we - kd choos ment put one whic is maduly restrictive \lore 
rototy pe a us To U . iso OLA cr¢ eTALL “7 ~ i square root ol | not further 
— + 5 cower ' : propo - Lity cle sf babl I} Is is i inlogo Is tO Lu complex 
onsta Ohm's Law experimentally umbers nad j 1 in quaternions The 
We we - IY Have \ irbith a : -_ s tol oul (; bbs Scull ind vector products ¢ dimensional 
SVvstel oltn irement units, but ilso uN oices ol implic tions 
space Hence \ 3p Cs ( ar ssigned arbitrary sO 
mes or dimensional symbols, and limensions WF. d yields [VW / 
ot the el ! spaces express products whil 
cluding reciprocals) of the indepé Ven. The T—d®F yields [7 } 
V ( 1 cle cj 0 a one itors ac) 
of the group of dimensions, and are also often calles T|=(W 
mish Basi ou out fers to be 
Piast uply ne Pos si 0 \ ind s ( | 1 1 
Purch ea dnne | i ( = 
e ( tity e Oo ippropriate to a give whit o Ss wit! our verb | deseription of torque as 
rol deduced by application of bas vork pel t angie 
Ws ( $ ol ition oF the basic equ ons Att Wt issivl / } cleterTm or the 
l mine dil , amma tes evrern rue su ement units, we still | bitrary 
nro mndiuct oO ‘ ors - CO} ho ( j jh Wi make these 
Oo leeb pplication of dimensional | €XP rv imphieitly b rl choice of {q] 
| or ti v quantity in an eq 0 
oa | Kor « mple, the force betwer wo current ele 
(‘hoos | | oO be imeric oft discard : : ae ~~ 
Wquely ( es Ss represt LION Spare I s s 
7 an i ae J meric. the converst iF "11 . r 
o lo true, and a dimensional label does no a 
uel tity rie aortq init \ Ko CX mpl I} irb choice ot | 1+ () defines the 
orqu cl ergy are usually comnu measureme unit of ecun erms of those of 
ain sion label, w t! vial attempt it distinetio fore dad dist ay Classical] | = ssigned to 
hy dis _ , bye rw “foot pounds (| woul bye | ‘ elec romagl tix ol , s ind 
f Tl s acthe liv distasteful. since torque ts . 
Worl ) | " or], The contusiol ses trol 
s> 7 erie lO ceonvenlencs or because 
\Maxy loc. « referred to “all quantities esset 
ti lly l sucha exponent ind exponent ils 








dimensions, yielding [J]=[m]'? [/]'? [4-'. But if we 
do not assign a space to T,, we can assign an arbi- 
trary space to J, yielding: '° 


r,] [F| [m] [2 
™ WP WPler 
This is the customary procedure associated with the 
MKSA system of units. There are four 
symbolic units (meter, kilogram, second, ampere 
associated with four basic dimensions, or generator 
spaces, but there are still only three basic measure- 
ment units. The measurement unit of current 
(ampere) is defined as a derived unit by assigning 


basic 


a measure to [,. This illustrates the confusion 
arising from the use of the word ‘unit’? without a 


modifier; the (measurement) unit of current is a 
derived unit, but the (symbolic) unit of current is a 
basic unit. Even worse is the substitution in the 
above of the name of a unit. The ampere (meas- 
urement unit) is derived; the ampere (symbolic unit 
is basic! 

Any proportionality constant appropriate to an 
experiment in vacuo can be thought of as a property 
of space. This interpretation is frequently em 
ploved when the phenomenon involves ‘‘action at a 
distance” and the proportionality constant depends 
upon the intervening medium. ‘Thus, in Ampere’s 
force law, the proportionality constant is customarily 
factored into two terms, uI,,, where u is a numeric 
(unity for vacuum) representing the differences 
among media. By virtue of a field interpretation of 
action at a distance, iv is thought of as a measure of 
the permeability of the medium to the field. Hence, 
we find the usage: 


uw is relative permeability 
ul, is “‘absolute”’ permeability 
l,, is the permeability of space. 


The term “permeability of space”’ implies a property 
of nature, experimentally determinable, but T,, is 
arbitrarily assignable! The measure, {T,,), defines 
the measurement unit of current. The only prop- 
erty of space involved is represented by the “hidden 
variable’, speed of light, which occurs in the con- 
straining relation between two otherwise arbitrary 
proportionality factors: TT, l. 

It is just as logical to interpret the constant, /?, in 
the ideal gas law, 

pV RT 

property of space. Elementary thermody- 
namics shows us that for any ideal gas, ?=¢ ( 
the difference between the specific heat at constant 
pressure and that at constant volume. Experi- 
mentally, PR is determined by a sequence of measure- 
ments of pV. z at lower and lower us cdlensities, 
finally extrapolated to zero density! Does this not 
yield (C,—C,) of a vacuum? 


as a 


If we con 


ider the 4" in the denominator to 


The constant R could have been arbitraril 
assigned thereby defining the measurement unit of 
temperature For reasons of convenience, howevel 


the unit of temperature was chosen as basic (det 
mined by a prototype standard . leay Ine fi is 
constant to be measured 


Rationalization 


Coulomb’s Law is not a eood Starting pot lol 
intrig 
The introduction of 42 into tl 
denominator is compensated by changes of T, and 1 
to keep force unaffected. 

We start with the 
tions relevant to the proportionatlities 


discussion of rationalization, despite the 


occurrence of 427° 


instead electromagnet t equa 


bu dra! 


f D-day 


In particular the magnetostatic field of a current 
loop can be expressed In terms ot a scalar pote ntial 


oalQ 


where Q 1s the solid angle subtended, at the pomt ol 


interest, by the current loop It is convenient to 
write 
q2 
© ke] 
Ss 
and 
VQ 
H—k/ 
\ 
where Ss represents the complete solid anel sul 


rounding a point, i.e., the solid angle subtended by 


any closed surface at any interior point 


The multivalued nature of the solid angle Q, alor 
anv contour linking the loop vields 
pva2-da=S 
and therefore 
pP H-dA=k/ 2 
We also have the geometric relatiol set ippe 


dix 2 


without prejudice as to the dimension of Q 


Hence 


unit or 
the angele subtended by thre loop Is 


*r- nda 


j 
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rationaliza 


is there is no 


Lit ! iy isure this relation finding 
£4) 
USS 


In bot svstems, eq ) requires insertion of the 
 f \ unit vector | | 
\ 


us au tactol it this quation is to be 
\ considered i quantit equatiol There iré three 
J CTLOLCeS 
rT C‘onsider both torms of eg 4 is measure 
1 
; equations ONL 
s ix b) Insert Ig explicitly in writing these equations 
| - 
: Require [To] | 
USstol Vy honrationalized equations require 
( now 


( ( makes [Q A i so that solid 
Ly ingle and plane angle have the same dimension and 
f \ | unit. 1 solid angle 

"2 


is Measured i 


ordinary radians 
Ol else solid radians and pla ( radius ine hot} 
j br hummer 
| The parallel analysis of 
‘S 
v-D=i 
™ v-Ddr=kg=QD- do 
\ Ilr 1 »* ‘ 
| I 
D | ev (-) 
S 
Th is ! }, fe Cuses solid nol Ss mensured nh stera 
adinns S Ir implies S br (st! na | sa unit vields hie same conclusions throughout and the 
ctor for accommodating dimensions results are compatible with the Maxwell equation 
\s po ( out b Per een thre most suitable ior Vv H 
terp oO ol | tr is | Lr =} S oO ry T) aim S10! al co clusio S 
I] | 
H A=SI 
4 I] Q| [J 
iz f | q 
lf w oul terpre S S 
\ could aqenuiy /] Wil HI] S ( re ol D Y 4 : 
solv i 1¢ t to CXpl t! eratlo Il Lo 
// I} Ss 1s Youne’s 1 rpret illol Ol thre effect 
? ol , ’ he « e] , t of soli ' ] bs 
oT ' . ot aie cae ¢ uni P solid angle | wis ot c ] B. Silsbe for many stimulating 
id would 1 solve | ! rn tas to Ww ethel q inh nd ve iabli dl scussions 
es oO wuhits ine ehal ect by rationatlizatiol ! 
vor of thre practi | ewpoiit Unfortu itely 
lo Ss viewpol we have already found if Notation 
‘ lr tor co ( | e 7 oO ize iorn 
~ s itico ~ ( re- 
O ; H A=/. is to b | , eis | o pl 
ql 0 ( !) omnes s i | 
| rive 
, I stract 
[VQ lV! : f 
H 7 
S tr (st! 
\ ' 
radi 
76 
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Appendix 1 
The meaning of (p)(s) cannot be directly inter- 
preted in terms of experimental conditions. 
Consider 
E=(1/2)mv? 


(E’)=(m)v)* 


The unit equation does not imply that a unit mass 
moving at unit velocity unit energy. 
Realizations of measurement units are inferred from 
measure equations. A system possesses or exhibits 
the measurement unit of energy, £,, when |) =1 


Now 


possesses 


BE m p}?/2 


so that /, is realized by a unit mass moving at a 


speed which is y2 times as big as v,, or by a Mass 
twice as big as m,, Moving at rv 
Appendix 2 
Let us choose the radian as the coherent unit of 


angle, so that the familiar measure equations 


6° 
. ) . —_ 
cos f ] ——' 11 

0 #7” 


sin @ 6- oe 
o. o. 


can be adopted, without change, as quantity equa 


tions. Then 
6.7)? 
cos 6 |— —~ + 
4 4 
sin 6 0:0 = { 
ede 


For each right hand side to represent an element in 
a single space, rather than in a mixture of 


we must have 


spaces, 
Ai l 


This implies that any variable z such that its space 
[x] satisfies [x]?=[1] is an allowable argument for the 


trigonometric functions. Thus, sin #: is also 
defined by the infinite series. 
Since (coherent = (radian radian )? | If we 
choose a noncoherent unit, Say, 
6 a radian 
then 
“ a 
and 
a ¢ i] ud 
radian)” ist I i 
elen dian’ 


This yields 


a a a 
sin @ ( @ : + ) 
ck 1). 
ad ‘ad 
Sihaé sinh @ , 
ie a 
so that 
Sih @;<(sin G7) SIN (6? raq/ Qa radian) sin | 6 
where (sin) is the unit vector in the space where sines 


are represented 

It follows that [sin] A|, space is 
unique unit For 
coherence, the symbolic unit of sine is the radian 


and if a 
vector, 


riven 


to have a sin i 


The consequences of Al > 1] ure mteresting Th 
vector area of a parallelogram is °r r and I 
volume of a parallelopiped is r r These 
vield 

Area /\7|@ 
Volume / 0 
which uppear odd but are compatible W 
CGrauss Divergence Theorem 
|v Vir=DV-do 

’ kes’ TI re 

nha Stokes hneorem 
| uxV-do=DV-ar 

Solid angie is oltten defined in terms of sp ! 

rea by 

re) iA 
or its equivalt I 

S2—sin 6 dé de 
With Area] A] and sin @| 0 0 | I 
formulations agree dimensionally and vield) (Q 

i radian The lack of distinction betwee! 

‘solid ingle” and “plane anel ’ is in ilovous ) 
lack ol distinetion betwee! plane ingle’ 
“numeric” when dé@—dl/r is used as a qu 
equation 

We can avoid all reference to area by consid 
a finite iit ute solid angele bounded by pli nes fg 
\t ene hy edge thie dihed il angle can bye expresst ‘| 
terms of its supplement, 7,, the “turning angle’’ of 


the normal to the edge It is readily shown (in ter 


of the « xeess of 
solid angle is 


by Which oF 


a spl erical triangl that 1 1 


ure of the proportional to t| e 


nu con plet turn 


BD hiskanaee U.P.T 


PT is the “uncompleted part ot one tun 
The general solid 
units to the corresponding plane angel 


convention for relating 
units is tl 


the measure of a tie Imnisph re in solid units is the same 
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This is an illustration of the fact that equations used 
in pl ysics need not be dimensio! illy homogeneous, 
when terms arise from the use of an artifice The 
form e” does not actually occur in physics; it is 
introduced with the convention that its real part 
alone or its Imaginary part alone represents a 
variable of interest Nonartificial « xponentials such 


us ¢ , have humeric exponents 
‘el r The turning-angle approac! to solid angle leads to 
' a formula for the solid angle subtended at a point by 


Hn irbitrary closed space curve 





RI where 7’ =dr/ds, and ds is any parameter of progres- 
SLO! round the curve 
8 Ul measure OF a cre n plane nits 0 
hemisphel 3 2n steradians or 360 solid decrees 


2 


‘} is for any unit of plane angel we cadefine the Appendix 2 


eresting to note tive sé itiol between 


2 | P 'T 1 ris wa “radians Althou | ole complete 
] 


‘his mopli turl Ss equivaient to 2r radians ( anv closed 
: | PT pian curve subtends 27 radians it an internal 
; ; po nt thre concept ol tur is topolog.\ il and not 
\) 1 o] . re metre rie rhtity ber OL Col plete Dusscs around any 
es oO oO closed pat! independent of our concep Ol angle 
: But is SOO! is we trv to mensure Iractiol al turns, 
PT \\ ist deni a metre concept tor the nterpolation 
But ¢ oa ; : ii ‘ nad this volves defining a 1 t lor ensuring the 
“ty — , rr op - | ; oe : poser orl ol i martial tur Wi Ist ilso speciiy thre 
ian a ial seiniins d J ls seinimans ocation o ¢ po taro a Wi I vi r onsidering 
| \ oO +1 Ox onl ole igo «lis lis! ri » rl I 
. , Co ( ire ul mit ro (i ong straight 
pelwe Olle Tier ql pla ! ’ , , ‘ / 
\ Ee a ; ¥3 | ela : eC, ~ the circuital law deriy on ot the magnetic 
s a held hent ong straight current ! ne conductor 
nea we a - “aug This I oby ous LLVOLN 2t Lillis If 
: ‘ nowe ( eonsietl 0 1) Osc ( rve We 
“her v9 ss ; 2 b ous ) { —_— ‘ et there are pomts on the wire for whucl 
| ) tri a I ) oO} ( t COS ~ fy} Jar ( ) ‘ i) er ~ ot rr illy 
Os Se iL ited W ( off rs ° 
du S10 lp dox. | owevel - olved ipp Pyle 
} ‘ , p ~— i \\ r) ras | electro iwnetic theors CO) éfe eneirctements 
oO . < Ke id enclosures unre roived wu thre basi 
eoncepts and equations Reso illo 0 total effects 
to parth effects associated wit! ylearne angles and 
, , solid a ries it ids to paradoxes isin the conventional 
oree law or current elements 
I aquatio piace no restrictiol oO ~ “he If turl is topological and not eUuri it should 
paradox comes from setting sin 2/2 | ther thar ot appear as a dimension or unit hus ““ampere- 
n 2/2—(sin 7 turn per meter” is a hvbrid ter “ampere pel 
iso tt st to note t il ¢ s not only metel should suffice The inclusion of “turn” n 
two-dimensional in the sense of having both real and this unit for #7 is also redundant, since linkage 
macmary parts but is also two-dimensional is relations and current paths cannot have fractional 
Compo! ni in two spaces turns in\ | kage olan pat Wit! x ere-Lurns 1s 
sini pi " nkauge Wit X aliperes 
4 eosd ) 
sf ] 
sin 6 A Paper 65B4—60 
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On the Range of a Fleet of Aircraft 
A. J. Goldman 
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for g<G,, where the caret denotes an omitted argu- 
ry . . . 
ment. The inductive hypothesis shows that 


M(y; Yi, Ya,» - + Ya-r) 


is a non-decreasing function of its last n—1 argu- 
ments over the set 
OSMmS#eS-.- Sy 
and so, using (2), (8) becomes 
M,(g) max (2+M(g—nRae; go, . . .. 9; (9 
rsg/nk 


for g<G,. The inductive hypothesis shows that the 
polygonal function 
P(z) r+M(g nr; Go, . « oo f LO 
of z, at all points in 0<g—nhr<—G,_, for which it 
has a derivative, has a negative derivative, and so 
P(x) is decreasing for 0< g—nRr< G, 
If g<G,_, then the maximum of P(x) required in 
(9) is therefore obtained by setting r=0, vielding 
M,,(g) M(G; G2, o « «9 JGuje (11 
But by the inductive hypothesis, the right side of 
(11) is given by the right side of (6) for g<@ 
Thus (6) is proved for g<G,_,. 
MBN ext suppose that G,,.<9<G@ 
function P(x) is decreasing for 0<g—nhr<G,,.,;, as 
noted above, and is increasing for g—nhr>G 
since by (7) and the inductive hypothesis the second 
summand in (10) is constant in this range. Thus 
the maximum of /(z) required in (9) is obtained by 


The polygonal 


setting g—nRr=G,_;, and (9) becomes 
fi G (j G 
M,(q)=” re Ge... : 
/ aR t 7 ad nk 
(G G G 
_ I + + are + 
Th m—1'1-2 n—2)(n—1 
R G4 G ios G 
n’'1-2 n—I1)n 
so that (6) is proved for G,.,<9g<G,. Finally, (7 


follows from (6) and the obvious fact that 


M,.(g) M,(G, ) if g G 
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3. Equal Fuel Capacities 


In this section it will be assumed that 
qd G for / 


and that the aircraft are so numbered that 


/ A, 7 >? i) 
It is convenient to define 
R 1/>> 7 for Lo 4 14 
The maximum range is denoted \7,(g), or mor 
fully by 
(G31 MC 


If q: nd, then 
which 


The maximum range is given by 


unique integel iy fo 


V if) G R R ° R R kG 
il ¢ rity Li 
We G(R,+R R fg >nG S 
In the optimal policy, if g<nG, the flight is bi 
with aircraft n, n—1 n—k If g>>nG, the 


nd vallons of fuel must be discarded and ihe 


( 
flight is begun with all » aireraft In any case 


the flight is begun with aircraft n, 3 g 


then the last-listed of these aircraft is abandone: 
when only tG gallons of fuel are left, the next-to-la 
is abandoned when only (¢ 1)G gallons remain, et« 
The proof of this result is just like that in 
preceding section, and is therefore omitted By wa) 
of orientation, it may be noted that the analog of 
9) is 
V7 qj nar J / qj R 


> ~« ’ 
Paper bobs ol 


Measurement of Wave Fronts Without a Reference 


Standard: 


Part 1. The Wave-Front-Shearing Interferometer 


James B. Saunders 


lune 26. 1961 


Introduction ol the separation of the points This irrangement 


and spacing of the reference points causes then 


vave-lront-shearmng iiterteromete! i, 2,0 laces, in the region ol overlappu ry wave fronts, to 
provides comparison Of a converging Wave front eomeme n pairs as is shown In firvure LA 
\ Sie ed —! Ol sell Su I Chie Phe separation oO} the reference points will be 
Wit ¢ ront-reversing mterleromete! I provides a defined as unitv for this discussio1 If the macni- 
COMLpPariso or either a converging or phihe Wave tude of the shear 1s unitv also. one adjustment ol 
ront Witt reversed uge OL itsell \ detailed the svstem (one set of fringes) is sufficient to eval 
deseription of at thod of analvzu the data Iron unte the deviations of the wave front at the chosen 
these mtertero eters, ¢ cept tor tle Spe i A ise Ol relere c¢ pots iro i Statist cally ehosen sphere 
revolutiol SViInimecrs has not bee published kor shear that exceeds unity more than one 
Consequent the absolute shape of a completely Iringe pattern 1.e., two ol ore shear values 
IKTOW Wiideé iIront co ld A be obt i! ed nitust bye ised Lo obtain nu SOLUTLLO 
Thi ethod of analvsis involves a mathematical In figure 1, P,(v=0, 1, 2 \V) represents the 
operatlol Chis operatiol is similas or the two \ | relerence points in one ol the images, W. 
lorementioned mterferometers lts deseriptior Is of the sheared wave front and /” the corresponding 
ore ClLOarly portraved its application to tle pomts im the othe lage WV’ Let ¢ see hig 1B 
wave-front-shearing interterometet WSI) than to equal to the deviation of the wave front at ?,, from 
‘ wave-front-reversing imterterometet! WRI thre corresponding pomt / on a reference circle, C, 
i oug the method evoived tron 1 study of the to be chosen late é equals the al orl between the 
atte There are. however, s ifficient differences i two unages, ¢ and C”. of the reference circle at their 
the two Operations to wari { separate treatments point of intersection and uw equals the distance from 
Lecording|\ part l whic follows will describe the mtersection point to P The distance, yu is 
the analysis of data from a WSI and part 2, to bi positive if /’) is below (in fig. 1B) the intersection 
published tel will describe the analvsis of dat of the two circles and negative if above The devia 
Irol 1 WRI Part 3 of this series will describ t101 is positive if 7, is on the concave side of the 
thie ibso te testing Of opt | flats This metho cirel and nevative if it lies on the convex side of 
~ ot new but has not been deseribed adequately the reference circle The distance from P, to P, (or 
n the literature from 7), to 7 measured along the circle, is always 
P positive and is represents | by 
art | The separation, S, (distance from 7, to 7 . 
l. Analysis of WSI Fringes . the two images of the peserenes circle at no 4 
pair of points, 7’, and 7’,.,, (assuming unit shear) 
‘he WSI will vield unique solutions of a wave | May be obtained from figure 2. The centers of the 
front. The assumption of svmmetry e. that it | two eircles are located at C, and ¢ The distance 
ive revolution svmmetr\ is unnecessary Any (,f equals 2h, P equals the radii of the two circles, 
onvergent wave front mav be tested along any p equals the distance from FE to 7 i. bisects the line 
cleaaum diameter of tlone the line of its intersection | CrC; and @ is the angle subtended by (y+ y) at £. 
with any plane that passes through the point of Applying the law of cosines (from trigonometry) to 
convergence Several equally spaced — referencs triangles AC’, 7 and £07, we obtain 
points are chosen along this line The line of ref 
erence pomts must bye parallel to the direction ol R / p ‘S Qh(y ‘ eos (90 7) and 
shear and the shear must be an integral multipl 


R f p Php cos (40 o 


239 











Ficure 1A Two images of a wave ont 


Figure 2 Relationship between parameters and othe 

On taking differences and solving for S, we obtain, 
S,=2h sin ¢. 

This equation may be applied rigorously, but if the 


focal ratio (radius of curvature of the circle divided 
by the diameter of the wave front) exceeds six, no 


Ts 
1. 


significant error is introduced by replacing sin @ b 
its approximate equivalent, (v+y)/?. Also, since t 
angle e will alwavs be quite small the value of 24 / 
approximately equal to R-e and the value of S 
eq | becomes 
. ) 

If we let ¢ represent the separation of the two 
wave fronts, Woand W’, at P,, it is seen from fis 
1B that 


J 


All quantities that represent optical distances a 
viven in units of the wavelength that is used sr 
cordingly, these quantities should be multiplies 
the wavelength to convert to standard units 
length The quantity (@, may be observed direct 
if white light is used to adjust the zero order to 


I 


a known point, 7, and then, using monochromati 
light of known wavelength observe the ordet! differ 
ence between points ?, and 7? It is inconvenient 
and unnecessary to use white light. Let Q, represent 
the unknown order of interference at an arbitrarily 
chosen point P., in the fringe pattern fig. 1A) and 
q the difference in order (number of fringes) betwee! 


this point and 2? The order of interference, Q,, at 
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I} ms mtroduces 


observation 


eretore equal to (, ( 


anothel unknown into the equations 0 


product 


but this mav be absorbed by repla ine the 
eu With a new paramete! such that 

€u () } 
With these values for (, and €u replaced ) the 
respect Ve equals eq becomes 


2 \ J 
epresents .\ 
ns namely \ | 0,'s 


equatlo Ss that contain 


additional 


two parameters € and Wi need. thre 
equations relating these \ unknowns Ia solu 
tio! s to be possibl It auppe ared that one could 
ise a second set of fringes formed with a shear of two 
its is now! i! i ie SA 77 Is would vield 
\ viditional equations of observations wit! 
0 vo new K1OW ~ two more parameters 
é BaL oO te wiced to thre noove-mentioned 
\ ; owns If \V is 6 or Jarger, we would 
nave is Ina eq mtions as u knowns ihe i sol ition 
wo I seen oO De DOSS byle Dut Ci be Show that 
ese equations are ot ¢ ely indepe ‘ 
| ( s (0) 0) nave sivtl iheanes | ( retel 
( cle f Show mourns iB is del ec \) 
circle Mma De ¢ ned b nres 0 1uionhs Ina il - 
} eo! | one I 1 rye ‘ ool ites OF ‘ 
“iso be ¢ ( is the o it passes throug et 
( Deo > ( ( cle ot retere ( Tri re full 
cle ec ar j Ons ol ond 0 Why 
thes, rire omlb the se o obs vatio 
equatio q (5 fron ~ e set o TES 
. {} ‘ t fy or mle solutiol 
0 ow 1! ilue of 2 Hlow ( 
est a parabol muro ol sta ( cif i 
ons represent the difference between the parabola 
j aspire compariso s to be Ihade betwee! 
he m sure deviations ol e wave tro FT| ‘ 
compute ( mations of i parabola trom -e rel 
the radius of cu iture of the circle must be OW 
wudius, /?, can be measured wlependent of 
tiie erterometet! ( deal nace point / 
re 2) or the mea Do tot converge ce ot e wave 
rol en Ie low il bhi its distances iron brh\ 
( On 1) t cal bye Tike iredl directt The eentel 
ol ‘ references circle | Te low med il ( ‘ 
mac 0 il ( elere ce emcl - 0 is ati 
eirel best fits the observe wave tro This 
ma Tye 0 most precisel 1 the metho ol leas 
squares ( iwlue of P? is the measured distances 
ron ( eal image po point of convergence) to 
i chose relterence point plus the computed Ccevia 
0 ol nis relerenes omt from the reference cirel 
| the test ol a concave murrol /? is approximately 
equal to the radius of the murrot n the test of a 
lens (simple or compound systen P is the distance 
irom thre pot ol convergence Image po oO the 
rae surface of the lens 
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2. Specifications for the Reference Circle 


There nr three relativel simple methods of 
defining the cirele of reference Kach method re- 
quires three equations ol conditio These con- 
ditions involve three well-known principles: namely, 
the method of coincidence, the me “i Ol averages, 
und the method of | ist squares 


2.1 


We ma define the retere rele shown 1n 
ives 1B and 2, bv requiring to pas rough 
} three ol the chosen reterence p Ss 7 Is, We 
might require hat it pass throug rhe two end 
po s P, and P plus some other point near the 
axis of the mirror (or lens If e points, 7 and 
PP? ur’ oO close to a dubbed « a bette fitting 
‘ cl Ss OF wned D\ requiring t to pass thro wh the 
vo points adjacent to 7’) ane Pp respectively 
Thes vo oices for reference ¢ le are made by 
equating e corresponding Oo zere Anv three 
ot the points can be chosen to cel e Ule circle and 
the results would be the deviations of the wave front 
it all chos reference points from the chosen circle 
For illustration we will choose the fi re pattern 
show! n figure 1A There are 8 reference points 
We will arbitrarily require the eference circle to 
pass | rougt points r P wna P This is done by 
requiring t it 
i) a) 
For this illustration eq (5) represent 7 equations 
ol obser itiol These ri equatio Ss. wie combined 
with thre condition equations t form a set of 10 
lineal equations Wit! 10 unknowns: and all s that 
ue not evaluated in eq (6) can be computed in 
terms of Known quantities the g,sa dy’s Table 1 
1 fperurs $ show thre results obtained by this 


2.2. The Method of Averages 


od of averages will usually vield a better 
; 


The 
fitting circle than the method o 
the 6,’s. This method 


\ reasonable correlation 


met 


nes bec nause 


averaging ol requires 


rnment of weights 


ot the 


an 
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TABLE 1. Deviations of a wave front from reference circles that 
are chosen as follows: Column 3, by the method of coincidence; 
column 4, by the method of averages; and column 5, by the 
method of least squares 


1 
| 
| 
| 


é, 
, de -_ 

M¢ MA LS 
0 0. 00 0.20 0.10 
l 2. 74 23 13 14 
2 2.10 15 ll 07 
3 1. 70 00 00 O7 
4 1.45 06 05 13 
5 1.30 06 04 Ol 
6 1.00 21 14 l 
7 0. 30 00 —. 16 Ot 
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VALUES 
icure 4 Graphical representation of the deviations of a wave 
front from reference circles that are chosen by the method of 
coincidence (MC), the method of averages (MA), and the 
method of least squares (LS). 


of image quality, formed by the measured wave front 
from an optical, image-forming system, with wave- 
front shape would normally require that the weights 
should be proportional to some function of the dis- 
tance from the reference points to the center of the 
circular aperture of the optical system. This would 
allow the assignment of equal weights to equal areas 
of the wave front. Equal weights will, however, be 
used in this discussion. 

In previous applications of the method of averages 
[5], the deviations of the observations from a curve, 
which is to be fitted to the observations by adjust- 
ment of parameters, are equated to functions of the 
parameters. These are then grouped into as many 
groups as there are parameters to be adjusted. The 
algebraic sum of the deviations in each group is 
equated to zero. When the functions of the pa- 
rameters (eq (5)) are substituted for the deviations, 
there are as Many equations as parameters and the 
parameters can be evaluated. Any manner of 
grouping the deviations will yield a solution. A 
very poor and sometimes useless solution will be 
obtained if the points that are allocated to each 
group are scattered among those of the other two 
groups. This author prefers to avoid any inter- 
mingling of the points of one group with those of other 
groups. 


The parameters of the reference circle shown in 
figure 1A are e, r and any one of the deviations, 

The parameter ¢ (see figs. 1B and 2), together with 
, determines the distance, A, from the center, ( 
of the reference circle to the point of convergence, £, 
of the wave front. The center of this circle is 
displaced from the ideal image point, # (mean point 
of convergence) for two reasons. One is that the 
axis of the mirror or lens cannot be adjusted to pass 
absolutely through the image of the source, and, 
secondly, the interference fringes would be too broad 
for reading when the wave front is nearly spherical 
Consequently, the angle « is adjusted to be large 
enough to provide fringes of a most convenient 
width for accurate readings 

The parameter, uw, together with 2, determines the 
direction of the point C,; from £, relative to the axis 
of the mirror. The parameter 7 is related to and 
replaces uw through eq (4). Consequently, the 
parameters 7 and ¢ define the center of the circle of 
reference relative to the wave front 

Any one of the 6,’s may be chosen as a parametet 
By adding 6, to the directly measurable distance 
from # to P, we have the radius 7? of the circle of 
reference. We will choose 6) for the third parameter 
because of simplification in the formulas that follow 

This relationship of the three parameters to the 
coordinates of the center of the circle and its radius 
is somewhat vague, involved, and perhaps not 
clearly explained, but it is sufficient to understand 
that they correspond to the three parameters that 
define the circle of reference 

If equal weights are to be assigned, there must be 
an equal number of points allocated to each of the 
three equations of condition. In the example chosen 
above (fig. 1A), the number of points are not inte 
grally divisible by three. However, we may considet 
that each of the observation equations represents 
three identical observations. The total number of 
observation equations is then 3.V and division by 3 
Is how possible. The three equations of condition 


are 
0 0 26 0 0 0 0 26 0 0 0) 7 
The set of equations represented by (5) may be 


replaced by the equivalent set 


where f, is a function of the three parameters 


e, and r. Any equation in set (8) is obtained by 
adding the first vy equations in set (5 The functions 
f,, may now be substituted for the corresponding 
6,’s in (7). It will be noted that there is no / 
corresponding to dp. This, however, does not present 


any difficulty because 4) is itself a parameter and 
may be left in eq 7) as is 

The three equations of (7), after replacing 6, by 
its equal, f/,, contain the three parameters as un 
knowns which may be evaluated. The computed 
values for the parameters are then substituted in 
eq (8) for evaluating the remaining 6,’s. The results 
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for the chosen example, shown in figure 1A, are 
given in table 1, column M.A, and are also represented 
graphically in figure 4 

The number of reference points, in the unit-shea 
determined by and dl- 
the element testing 


the shear angle 
tested For 
during the polishing 
the optician is only inter 
in departures from a figure. The 
easiest method of computing this is, therefore, the 
most practical. The shear may be adjusted to such 
a value that the resultant number of points is some 
multiple of three. The computation is a minimum 
when the number of points is the same in each group 


method, is 
ameter of being 
elements or lenses 


figuring operations 


mirrors 
and 


ested desired 


2.3. The Method of Least Squares 
The 


squares 


method of evaluating a wave front bv least 


is more ¢ laborate than the others but would 
desired 


whet highest 


method, the 


usually be used 


In applyvi Y 


pre¢ ision is 
! 


this sum ot the squared 


deviations ( >> 6 ) is to be minimized with respect 


to the three parameters required to define the 
reference sphere 

This application of least squares differs from the 
observations 
and the number 
number of 


the é’s 


usual method of least squares in the Lhe 


ure assumed to hye free from error 


exceed thre 
namely 
obtained even if 


of observations needs not 


unknowns whose values are sought 
In fact, a unique 


the number of 


solution mav be 


; 


observations Is less than the 


number of ( hosen reference 


one 
pomts 


To minimize the sum of the squared dev latvLIOnS of 


the wave tront from the reference sphere we require 
that 

\ \ 

> / ee. pe mini i! +e] 


Tl is is effected by equating to zero the differentials 


ol eq UY) will respect Lo é and do, respectively On 
performing the differentiation and dividing throug! 
bv 2. we obtar 

\ f v >i v >i 

a eS >: ¢, Y=0 10 

Kad ¥ Knead Pe — ay 


The pal tial differentials, as obtained from eq (S), are 


VO} , > 
| 0.1.2 \ | 
06 
f 11 
)/ , OF 4 . 
) Bate ye y 2 A 
F 


se differentials into (10 we have 
reference points fig. 1A 


On substituting the 
lor the example of 8 


Qn combining sets of eqs (8) and (12), we again 
have 10 equations with 10 Using the 
same set of fringes fig ] \ a8 Was above for 
illustration, we obtain for the 6, ’s the values given 
in column LS of table 1 These 
plotted, along with those obtained by the other two 


in figure 4 


unknowns 
used 
also 


V alues ure 


] 


methods 


3. Sensitivity 


The sensitivity of the WSI varies with the absolute 
If the shear is maintained 
the sensitivity of the observations increases 


the number of 


magnitude of the shear. 
at unity 
points because the 
The reduction of 
the number of points is, of course, limited 
In testing surfaces that depart very greatly from 
a sphere, as for a large aperture and large aperture 


with a decrease in 
unit, as defined above, increases 


desirable Lo 
such a murror 


the W SI. 


to-focus-ratio parabolic murror, it is 

sensitivity. Otherwise 
tested interferometrically 
number of points is increased, the absolute 
smaller the 
shear relative to aperture, the smaller can the optical 
path difference be made By reducing the absolute 
magnitude of the shear the maximum optical path 
difference may be reduced at will and, consequently 


reduce the 
cannot be will 
As the 


spacing of the points decreases The 


the fringes may be adjusted to any desired width 
The sensitivity may also be increased, without 

reducing the number of reference points, by using 

two or more shear values that are greater than 


unity hig > This may vield more equations than 


there are unknowns but averaging of computed 
points Is alwavs acceptable Each additional set of 
Iringes obtained with a different shear value, 


introduces two additional unknowns—new values 


for e and ? 


4. Symmetrical Surfaces 


venerated mechanically 
renerally have about the 
optical axis. If this is known and the axis is known 
to be at the center of the mirror (assuming circular 
apertures the known to be 
svmmetrical because the center of the reference circle 
Two other con- 
the 
symmetry 1s 


Optical surfaces that are 


revolution symmetry 


then deviations are 
assumed to lie on the axis 

then sufficient to 
revolution 


can be 


ditions are define reference 


The 
defined by the spec ification 


circl condition for 


ri 6 13) 
and the condition equations 6b 7) and 12 be- 
come respectively, 

0 0) 14 

} SY 2 6 2 0 15) 
rie 
ra) 0 6 {) 
; ahs. . 16 
216 156,-+-116.+-9 0 
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The third equation in each of the sets (6), (7), and 5. References 
(11) becomes identical, respectively, to the first | 5.) wey Bates. A Wavefront Shearing Interferometer. Pro 
equation of the set. Phys. Soc. 59, 940 (1947 


The assumption of symmetry (eq (13)) reduces | [?! BR. L. Drew, A Simplified Shearing Interferometer, Pro 
ri - Phys. Soc. (London) sec. B, 64, 1005 (1951 


the number of 4,’s by approximately 'sN (exactly }sN | 13) Andre Danjon, Mach Interferometer for the Measurem 

if N is an even number) and reduces the number of of Fringe Visibility; Arthur Beers, Vistas in Astronot 

condition equations by one. If N is larger than 2, | |, ,! Ned ee eee is Roel Mili 
the number of equations exceeds the number of un- ing Interferometer, p. 227 (Pergamon Press, 1960 
knowns, for a unit shear arrangement. Conse- | [5] J. B. Searborough, Numerical Mathematical Analysis, | 
quently, larger shears (2 or more units) may be used wee (The cohns Hephine Prem, Baltimore, Md. (1000 
with increased sensitivity, as explained above Paper 65B4-62 
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On the Evaluation of the Function 


\ ( ‘ay 
) 


for Real Values of i’ 


Wolfgang Borsch-Supan 


lune 26. 106] 


a. Introduction ) [ ne We obt ill 


ordet! oO dese! t! P - Pe 
Slo Lio 0 Spectt! adisti Dutio Ol ¢ J . 
H lese | | 
i ves Shor ‘ Scrip ) 0 ~ eto oO! CO | (1 j ) | | i )) 
) 
Ss] ! ~ r ld L plo re ~ plot = 
) I } ol = ( ~ eul \ : . — 
eg aes ‘or large positive values of \ this vertical line i 
ePSper oO 
Por ti ake the complex plane passes near the singularity of the 
( j eal e 1) ition Oo e 4 Oss . 
, } . ‘ , ntegvgrand at the orivt | ] ermore Live pl ist 
0 Cel ) of particles Cel edu ; rr 
. ‘ | } ‘ - ¢ . ‘ j angie ot the integrand varies el last Therefore 
( sto be ipu dand orappro ed , 
; : this pat of integration 1s not table lor umerical 
) i ae | i ~ ree i*- ‘> { 
’ = ile The eorrespondr erTrol 4 lidis oO! — - 3 . , Bt SE ¢ - cecreasing | 
\ th : the left half-plane for these values of \, we pull the 
v) \ les 0 i. ithe ore ( ci A 
ath of integration from to the left t 
, , P n fn , n {j ‘ _ a ev LbLiO we ( i = , 
" ne ‘ +} But the sir larity at the origin prevents a further 
rie sfO! 0 ( eC op . 
re movement of the middle part of the path of egra 
Oo 
| : P tio! s ce we have a bran Do tL a the origin 
‘ lifts re ¢ ‘ ‘ ollo 1) iwrauphs . 
we mut a branch cut atong the hegatis real axis 


‘ ad transtorm the integ! linto a ntegral along tl Ss 
2. Transformation of Path of Integration v. 


prane cut iron tougo | ( OW eC! 


As pointed out by Lan valuat 0 to on the upper side. Thy 


| vrand ditters 
sn +, een iin al gta b Dv tl difference of 27 1n the imaginal part of In w 
choo ’ hy ly p ‘ ly 7, | SO We 
1) oO oO 
—- - The sing ilarityv at the origil : win iltiplied bv a 
) sath factor of the ordel il compare with the original 
| integrand Furthermors thre tegrand decreases 
the faster the larger X is 
Phi ! hy ral 0 ed oO real one It turns out that there is a i inge in which 
I) mtrodue | / 3 a ew ndependent both formulas are of practical use but formula (2 
Variable and by combining the co plex conjugate is preli rable for large negative A, and formula (3) for 
ilues of the egrand for co rl x¥ conpugate Values large positive \ At the same time these formulae 


provide asymptotic formulas for \), given already 
by Landau, but without anv remainder terms 
eTro!l bounds We extend Landa is co 


in order to get such error bounds 
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siderations 








3. Asymptotic Formula for Large Negative } | 


For large negative \ the asymptotic behavior of 
¢(A) is obtained by applying the saddle-point method 
to eq (2). The integrand can be written in the form 


exp {os(¥)} - { og (¥)}. 


tO 


where 
4 In (i+#@)—+t arctan ¢, 


g(t)=t(4 In (1+#)—1)+arctan ¢. 

The following inequalities can be proved by elemen- 
tary considerations which are too lengthy to be repro- 
duced here in detail: 


fort >0. te 


, . t? 
Replacing f(t) by —: 


> 


and g(t) by 0 vields the asymp- 
totic expression 


o(A) ~ . exp ( a )dy — 
7 Jo oO \2r 
which was given by Landau. 
Let us now look for a lower bound of #(A) or for a 
lower bound of the integrand. We replace f(t) by 
f?/2, cos og by 1—(oeg)?/2 and then g(t) by #/6 
Integration \ ields the lower bound 


-( ' Sie } (9) 


[t is a little more complicated to get an upper 
bound, since different inequalities are to be used for 
small ¢ and for large ¢. We replace cos og by 1 for 


(nr) >>7 
\2r 


all ¢, f(t) by the upper bounds given by (4a) for 
t<t and given by (4b) for t>t, respectively. The 


quantity f) is chosen later, depending upon ¢. For 
the sake of closed integration, the expression exp 
4/12} is replaced by the upper bound 1+ (exp 
ti/12 1)tt for t<t). Now the integration can 
be carried out, provided to fy< 6, and we obtain 


en? (‘so 2 = y' 
rA)< ‘X —- 1+je7o**—1| = 
avs =f r{ az ; oti 
tf, # 
a ere 
t 


ty 


> 12 


Since the integrand of the remaining integral is 
positive for y > ot, too, we replace of by and obtain 


The value of f) should be chosen in such a way that 
the sum of the two last minimum. But 
the equation determining this minimum for general 

So we restrict ourselves to a 


terms 1S a 


is rather complicated 
rather crude discussion and to the 
The second term in parentheses is minimized for a 
fixed o by letting ofi—0 Then this term 
1/(4o But the third term infinity 

Now let us put of, equal to some constant, sav |e 
Then our inequality has the form 


cause of large 


OCS to 
then 


eroes to 


for convenience. 


valid for every \ and for every positive e< 30 \p 
parently, for o the last term goes much faster to 
zero than the middle term 
remains large 


So one can take e small 


as long as o-e Then the first term 


approaches 1/(4e 


4. Asymptotic Formula for Large Positive \ 


The asvimptotic behavior of @(X) for large: positive 
h can be seen from (3 In a manner similar to 
Landau’s, we transform the independent variable \ 


and the integration variable v. The reason is that 
the interval where the integrand is essentially differ 
ent from zero and contributes to the integral essen 


tially should be independent of \ for the new integra 


tion variable ¢. Such a transformation is given by 
A=w+lIn w+B 
f 
wr 


where the new variables are w and ¢, and B is a 
constant chosen below SO We get 
] of 
oir é'¢ resin f 
Tw, uw 


Tavlor expansion of the second exponential factor 
vields 


l rt l ae 
d(A é Sin dt ( t(In ¢ B) sin / 
w, u ru. u 
l : . . »— wt in fre 
+ e~‘t?(In t 2)? sin é if 
2ru w 


where 0< 6=0(t.w)< 1. 
In order to obtain a remainder term of the 
‘ My 

Ww we would like to choose B in such a way 


ordet 


that 
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ie same goal can be reached 


depend upon wu 
develops the 


Tavlor series too 


that the part of the second integral cor espondin 
Taylor expansion of sir 


the first term 1 
ermined by thre equ 


especially eq 


i similar 
Lie 
value 


Wi replace Live 


al can 


lv we have thi asvin ptoti formu 


to have the right hand sides 
Tl erelore one needs 


terms of A The 


It is more convenient 
9) in terms of A 


ot (S) and 
easily con putable estimates ol 
uncel of w due to this estimation 
uncertainty in the right hand side of (8), which in 
expressed bv a corresponding change of 
8 Therefore, in addition, 


\ Using the 


induces an 


rtaint\ 
turn can b 
Ain the left hand side of 

i an estimate of d’ 


we need 


value 


Ind = 10 


as an approximation for w, we obtain the inequalities 
I j 
u T 
ln X B exp 0.926/(u ly B 
p In wu Bw \ 
- 
, Lf? (u l 
| are valid provided 2.02 
5. Comments on the Table 
The numerical evaluation of Z mi ) lor 
given value of A requires the introduction of finite 
step lengths and of finite integratlo limits 
far as finite integration limits are concerned 
errors introduced by this approximat on can be esti- 
mate rivorousl\ The errors due tothe use of finite 
step engths and Sunpsons rule mre mucti more 
nvolved if expressed in terms ol easil computable 
is qua lies Repeated estimates would lead 
pool inal estimates Therefore, a mixed theoretic: 
show experimental way of guessil these errors 
In adopted These vuesses and some test runs 
alved step lengths allow e assumption that 
relative error of the values of @(\) as given in 
table below do not exceed the amounts shown below 
Together with the estimates of errors due 
ntegration limits, one obtains 
and 
LS LO il | AS2 
elat n table 5.210 f 2<A 
7 1O { () \ 
s 
For values of \- t the formulas (5 6), and 
vive the value of @(A) with a relative error 
ceeding 1.5 percent For values of X 100 the for- 
a mulas (11 10), and (9) provide values of @(A 
a relative error not exceeding 0.02 percent 
is almost 


The notation of the numbers 
yr. 1.25 f 04 means 1.25 LO 
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Table of the Function d(Xr 


aN (A) r (A \ or \ o(r 
4. 000 3. 3900 F-—09 2. 500 9. 6369E-03 1. 000 1 5140E-01 0. 500 1. 6524E_01 
3. 975 5. 4974E-09 2. 475 1. 0633 E-02 0 975 1. 5335E-01 0. 525 1 §6433E-01 
3. 950 8. 8066E-09 2. 450 1. 1700E-02 0. 950 1. 5523E-01 0. 550 1. 6341E-01 
3. 925 1. 3940E—-08 2. 425 1. 2840E-02 0. 925 l. 5705E-0 0. 575 1. 6247E-O01 
3. 900 2. ISLIE-O8 2. 400 1. 4054E-02 0. 900 1. 5880E-01 0. 600 1. 6152F_01 
3. 875 3. 3740E-08 2. 375 1. 5344E-02 0. 875 |. 6048E_01 0. 625 1. 6056E-01 
3. 850 5. 1619E-08 2. 350 1. 6712E-02 0. 850 1. 6208E-01 0 650 1. 5959E_ 0] 
3. 825 7. 8123E-—08 2 325 1. S157E-02 0. 825 1. 6362E-01 0. 675 1. 5861 E-01 
3. 800 1. 1700 E-07 2 300 1. 9682 E-02 0. S00 1. 65090F_01 0. 700 1 5762E_01 
3. 775 1. 7342E-—07 2. 275 2. 1285E-02 0. 775 1. HH49E-0O1 0. 725 1. 5661 E-0 
3. 750 2 5449EF-07 2 250 2 29068E--02 0. 750 1 6781E_01 0. 750 1. 5561F-0 
3. 725 3. 6984E-07 2. 225 2. 4731 E--02 0. 725 1. 6907E-01 0. 775 |. 5459EF-0 
3. 700 5. 3235E-—07 2. 200 2 6571E-02 0. 700 1. 7026EF-01 0 800 1. 53457E_0 
3. 675 7. 5919E-07 2.175 2. S490E-02 0. 675 1. 7138E_O1 0. 825 5254 0 
3. 650 1. O729E—06 2 150 3. 0486E-02 0. 650 1. 7243E-01 0. 850 1. 5150} 
3. 625 1. 5029E--06 2. 125 5. 2558 E-02 0. 625 1. 7341E-01 O. 875 1. 5046E-( 
3. 60 2. OS72E-—06 2 100 3. 4703 E02 0. 600 1. 7433E-01 0. 900 1. 4942] 
3. 575 2 &§743E-06 2 O75 § 6921-02 0. 575 1. 7T518F_O1 0. O95 1. 4837E-01 
3. 550 3. 9259-06 2 050 § 9200-02? 0. 550 1 7T596E_O} 0 950 1732] 
3. 525 5. 3195E-06 2. 025 1. 1565E-02 0. 525 1. 7668E-01 0. 975 1. 4627E-01 
3. 500 7. I51SE-06 2. 000 1. 3985E-02 0. 500 l. 7734E-01 1 ooo 15211 
3. 475 9. 5423 E-06 1 975 1. 64681 02 0. 475 1. 7793E-O01 10) | 1097] { 
3. 450 1. 2638E-—05 1. 950 1. 9010E-02 0. 450 l. 7T846E-0 1 200 1. 36731 
3. 425 1. 6616E-05 1. 925 5. 1608E-02 0 425 1. 7T893E_-01 ie 
3, 400 2. 1694E-05 1. 900 5. 4259E-02 0. 400 | 7O34E-01 mn rE 
. 375 2. 8128E-05 1. 875 5. 6958E-02 0. 375 1. 7970E-01 oo =oee-0 
3. 350 3. 6227E-05 1. 850 5. 9703E-02 0 350 7999E- O01 L. OOO ane + 
3 325 1 6353E-05 1. 825 6. 2488E-02 0. 325 |. 8023E-01 > 1. 2018E-0 
3. 300 5. 8931E-05 1. 800 6. 5311 E-02 0. 300 1 SO42E-0 1. (UU 1. 1622E-0 
3. 275 7. 4458E-05 l. 775 6. SI67E-02 0. 275 1. SO55E-01 
S00 l. 1235E 

3. 250 9. 3505E-05 1. 750 7. 1052E-02 0. 250 1 8068E_01 1. 900 1. OSSSE- 0 
3. 225 1. 1673E-04 1. 725 7. 3961 E-02 0. 225 1 S066E-( ». OOO 1. 0492E-0 
3. 200 |. 4489E-04 1. 700 7. 6891 E02 0. 200 lL. S064E-0] 2? 100 1. O1236E-0 
3.175 1. 7T882E-04 1. 675 7. 9837E-02 0. 175 1. 8057E_01 2. 200 9. 7911 E-02 
3. 150 2. 1949E-04 1. 650 8. 2795E-02 0. 150 1. S046E_0 

2? 300 » 45731 2 
3. 125 2. 6797 E-04 1. 625 8. 5760E-02 0 195 1 S030EF-01 > 400 9 1346-09 
3. 100 3. 2545E-04 1. 600 &. §720F 02 0 100 1 S0090F- 0] 2? 500 & S230E-02 
3. O75 3. 9326E-04 1. 575 9. 1696E-02 0. O75 1. 7O84E-01 2. 600 8, 5223E-—0 
3. 050 4. 7282E-—04 1. 550 9, 4658E-02 0. O50 1 7O56E-01 2 700 & 2394-09 
3. 025 5. 6574E-04 1. 525 9. 7611E-02 0. 025 1. 7923E-01 

2? S00 7. 9531 E-02 
3. OOO 6. 7372E-04 1. 500 1. OO55E-01 0. 000 TSS6E_O1 > O00 7 6842E-0 
2. 975 7. D861 E-04 1. 475 1. 0O347E-01 0. 025 7TR45E-O1 2 O00 7 4953-0 
2. 950 9. 4241E-04 1. 450 1. 0637E-01 0. 050 1. 7TSOLE-O} 2 100 7 1763E_-0 
2. 925 1. 1072E-03 1. 425 1. 0925E-01 0. O75 1. 7754E-O1 ; 200 6. O368E_0 
2. 900 1. 2953E-03 1. 400 1. 1210E-01 0. 100 1. 7703E-01 
2. 875 1. 5090E-03 1. 375 1. 1491E-01 0. 125 |. 7648E-01 oo 6. 7066E-02 
2. 850 1. 7509E-03 1. 350 1. 1769E-01 0. 150 |. 7591 E-01 a 6. 4853E-0 
2. 825 2. 0234E-03 1. 325 1. 2043E-01 0. 175 | 7530E-01 . 500 6. 27251 02 
2. 800 2. 3293E-03 1. 300 1. 2313E-01 0. 200 1. 7467E-01 3. OUU 6. O6STE-02 
2.775 2. 6714E-03 1. 275 |. 2579E-01 0. 225 |. 7401E-01 +7 5. 8717E-02 
2. 750 3. 0524E-03 1. 250 1. 2840E-01 0. 250 | 7332E-01 3. 800 5. 6830E-02 
2. 725 3. 4752E-03 1. 225 1. 3096E-01 0. 275 1. 7261E_0] ». GOO 5. 5016-02 
2. 700 3. 9428E-03 1. 200 1. 3346E_01 0. 300 1. TISSE_O1 +. OOO ». 3274-02 
2. 675 1. 4581 E-03 1. 175 |. 3591-01 0. 325 1. 7112E-01 i. 100 5. 1599E-02 
2. 650 5. 0239E-03 1. 150 1. 3831 E-01 0. 350 |. 7033E-01 4. 200 1. 9990E-02 
2. 625 5. 6431 E-03 1. 125 | 4065E-01 0. 375 1. 6953E-01 1. 300 t. 8444E-02 
2. 600 6. 3186E-03 1. 100 1. 42902E-01 0. 400 1. 6871 E-01 1 400 t 6959F-02 
2. 575 7. 0532E-—03 1. O75 1. 4514E-01 0. 425 1. O787E-O1 t 500 t 5530E-02 
2. 550 7. 8494E-03 1. 050 1. 4729E-01 0. 450 1. 6701E-01 t. 600 1 4156E-02 
2. 525 8. 7OORE-03 1. O25 1. 4937E-01 0. 475 1. 6613E-01 t 700 1. 2836E-02 
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Table of the Function o(r Continue 


X (A) nN (A » o( » or 
37. 200 8. 3033E-04 45. 600 5. 4406E-04 59. 000 3. 1903E-04 80. 000 1. 7017 E-04 
37. 400 8. 2113E-04 45. 800 5. 3914E-04 59. 500 3. 1351E-04 80. 500 1. 6800E—-04 
37. 600 & 1208E-04 46. 000 5. 3429E-04 60. 000 3. OS13E-04 81. 000 1. 6588E-04 
37. 800 8. O318SE-04 46. 200 5. 2951 E-04 60. 500 3. O289E-04 81. 500 1. 6379E-04 
38. 000 7. 9442E-04 46. 400 5. 2478E-04 61. 000 2. 9778E-04 82. 000 1. 6174E—-04 
38. 200 7. 8581 E-04 46. 600 5. 2012E-04 61. 500 2. 9280E-04 82. 500 1. 5973E-—04 
38. 400 7. 7733E-04 16. 800 5. 1552E-04 62. 000 2. 8794E-04 83. 000 1. 5775E-04 
38. 600 7. 6898E-04 17. 000 5. 1098E-04 62. 500 2. 8320E-04 83. 500 1. 5581 E-04 
38. 800 7. 6077E-04 17. 200 5. 0650E-04 63. 000 2. 7857E-—04 84. 000 1. 5391 E-04 
39. 000 7. 5269E-—04 47. 400 5. 0208E-04 63. 500 2. 7406E-04 84. 500 1. 5204E—04 
39. 200 7. 4474E-04 — — 64. 000 2. 6965E-04 85. 000 1. 5021E-04 
39. 400 7. 3691 E-04 ~ noe oki a 64. 500 2. 6535E-04 85. 500 1. 4841-04 
39. 600 7. 2920E-04 18 000 1 2915E_04 65. 000 2. 6115E-04 86. 000 1. 4664E-04 
39. 800 7. 2161E-04 418 200 1 84905E-04 65. 500 2. 5705E-04 86. 500 1. 4490E-04 
40. 000 7. 1413E-04 18 400 + 8080E-04 66. 000 2. 5304E-04 87. 000 1. 4319E-04 
40. 200 7. 0677E-04 ; Tan 66. 500 2. 4913E-04 87. 500 1. 4151 B-04 
40. 400 6. 9953E-04 48. 600 4. 7671-04 67. 000 2. 4531E-04 88. 000 1. 3987E-04 
40. 600 6. 9239-04 $8. 800 1. 7267E-04 67. 500 2. 4157E-04 88. 500 1. 3825E-04 
40. 800 6. 8536E-04 49. OOO 4. 6868E-04 68. 000 2. 3792E-04 89. 000 1. 3665E-04 
41. 000 6. 7844E-04 19. 200 1. 6474-04 68. 500 2. 3435E-04 89. 500 1. 3509E-04 

19. 400 1. 6084E-04 

41. 200 6. 7161E-04 69. 000 2. 3085E-04 90. 000 1. 3355E-04 
$1. 400 6. 6489E-04 49. 600 $ 5700E-04 69. 500 2 2744E-04 90. 500 1. 3204E-04 
41. 600 6. 5827E—04 49. 800 4. 5320E-04 70. 000 2 2410E-04 91. 000 1. 3055E-04 
41. 800 6. 5175E-04 50. 000 1. 4945E-04 70. 500 2 2083E-04 91. 500 1. 2009F-04 
42. 000 6. 4532E-04 71. 000 2 1764E-04 92. 000 1. 2765E-04 

50. 500 1 4028E-04 

42. 200 6 3899E 04 51. 000 $ 3138E-04 71. 500 2 1451E-04 92. 500 1. 2624E-04 
$2. 400 6. 3275E O4 51. 500 $ 2975E-04 72. 000 2? 1145E-04 93. 000 1. 2485E-04 
42. 600 6 2660E 04 52. 000 t 1437E-04 72. 500 2? OS45E-04 93. 500 lL. 2348E-04 
$2. 800 6. 2053E-04 52. 500 1 0623E-04 73. 000 2. 0552E-04 94. 000 1. 2213E-04 
43. 000 6. 1456E-04 73. 500 2? 0264E-04 94. 500 1. 2081E-04 
£3. 200 6. O867E-—-04 53 000 3 eer a 04 74. 000 1. 9VORS3E 04 95. 000 1. 1951 E-04 
43. 400 6. 0286E-04 sy ne : ene os 74. 500 1 9708E-04 95. 500 1. 1823E-04 
‘ . - "19R o < 5 4 Ji =e P . < 606 OA 
13. 600 5 97 13E 04 54. 500 3 7596E-04 75 000 | 438K 04 #5. OOO | 169 K-04 
43. 800 5. 9149E-04 Re On) 3 6202K-04 75. 500 1. 9173E-04 o6, 500 1. 1572E-04 
14. 000 5. 8592E-04 _ _ 76. 000 1. 8914E-04 97. 000 1. 1450E-04 
44. 200 5. 8043F-04 55 500 3 6207E 04 76. 500 lL. 8660E-04 97. 500 1. 1330E-04 
44. 400 5. 7502E-04 36. OOO 3. S541 E04 77. 000 1. S411E-04 98. 000 1. 1212E-04 
44. 600 5. 6968E—04 ob 500 3 I893E 04 77. 500 1. SIGSE-04 98. 500 1. 1095E-04 
44. 800 5. 6441 E-04 a7 000 3 12631 O4 78. 000 1. 7T928E-04 99. 000 1. OOROE-04 
15. 000 5. 5922E-04 57. 500 3. 3649E-04 78. 500 1. 7694E-04 99. 500 1. O867E-04 
45. 200 5. 5410E-04 58. 000 3. 3052E-04 79. 000 1. 7464E-04 100. 000 1. O756E-04 
45. 400 5. 4905 E-04 58. 500 3. 2470E-04 79. 500 1. 7239E-04 

The numerical computations were done on the 6. References 
National Bureau of Standards’ 704 computer. The 
author acknowledges the help of Miss R. Zucker, | [1] L. Lar de ; On the energy loss of fast particles 
7 . . a saan ois 8 tion, . *hvsies 8 (1944), 201-205 
who programed the computations. 19] A. Erdelyi, W. Magnus, F’ Oberhettineer, and F. G. 7 
comi, Higher transcendental functions, Vol. I, Ne 
York-Toronto London. 1953 
= 3] k B Hildeb and Introduction to numeri il tal’ 
(Paper 65B4 63) New York-Toronto-Londor LO56 
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Analyticity and Probability Properties of 
One-Dimensional Brownian Motion 
Abolghassem Ghaffari 
Ju 11, 196] 

Phe m il dime il Brov notio govel | 
Chapman- he ogoroff funct ' lation for! f serie of produ f ‘ 

Laguerre pol e described. Some at il and probabilistie propert 
, liseussed and ire different from t 
1 N. BK y fT r} I rf bat functi ‘ 1 
\ derivat respond ' tiff quation 
1. Intreducticn 
Che mathematical study of Brownian motior ind diffusion phenomena tive se Ol a 
oO form fluid led S. Chapman [1 to introduce the functional equation 
j ; . ; j j 
B. Hostinsky [11] in studvinge the metional equation of \L. Smoluchows! 19] for the 
Sure problem observed that the unctional equation oO ( hapman Was a gene! il on of the 
unectional equation ol Smoluchowski and found the reneral solution of the Ch pmal equator 
| ising the method of integration bv substitutions due to V. Volterra 
The same author [11] puts the functional equation of Chapman in the more symmetric 
orn 
I fy, r,t, f t, tedda 2 
which is equivalent to that of Chapman for the case where the space in whiel e Brownial 
motion ol diffus nm phenomena take place has only one dimension 
Qn the other hand A. N Kolmogorofl 17] has considered directly a very reneral fune 
ional equation valid for domain of any number of dimensions, which reduces to eq (2) in 
e case of a single dimension and which can be considered as a very general form of one of 
the conditions which occurs in the problem of chain probabilities studied by Markoff 
Moreover, A. Kolmogoroff has given particular solutions of eq (2) and has shown that 
( wide classes of his solutions satisfy linear partial differential equations of the second 
order and pal tholie ty pe and rene! ilized thre hunetior il equation ol hapman Irom the spec ial 
problems of Brownian motion and diffusion pl enome! to cham probabilities The eq (2) 1s 
commonly ealled the Cl upma Koln ogoroll functional equator 
Nevertheless, the most general solution of the Chapman-Kolmogoroff equation has not 
bre CLV ML. Fréchet, however, has given very general solutions having forms different 
om the forms give by the preceding authors 
M.. Fréchet has found two m ods: the first [6] is based on a procedure used mathe 
i" eal physics d the second 7| ol reducing the study of the most eneral solutions which 
ire doubly square ntegrabl to the corresponding study in the discontinuous cas The two 
methods provide two distu et torms whose equivalence is either obvious nor cert 
The most general solutions of the Chapman-Kolmogoroff equation in one-dimensional 
pace Tay heer Clvel by \I Fréchet id the present ithor In Iact the ease ol the finite 
iterval O27 has been treate | hy \I Fréchet t) il dl the case ot a region infinite both 
directions has been investigated in detail [S, 10 








In this paper some analytical and probabilistic properties of the most general solutions 
of one-dimensional Brownian motion in the interval (0,0) are described. The corresponding 
parabolic partial differential equation, which constitutes the fundamental equat’on of Brownian 
motion and diffusion theory, is obtained, and the limiting behavior of the solutions at infinity 


is analyzed. 
2. Statement of the Problem and Formulation of the Fundamental Results 


Let us consider a free spherical Brownian particle (i.e., in the absence of an external field 
force) depending on a finite number » of parameters and which can assume all states (, forming 
a certain set V (which one can consider as a region of »-dimensional Euclidean space). One 
assumes then that there exists a probability w(.M,s; vf) that the Brownian particle passes 
from state VM at the instant s to any one of the states of a set vy at the later instant ¢ 

The probability function w(.M,s; »,f) is, by virtue of the total probability theorem, an 


additive function of the set v. It is clear that 


w(.M,s; vs) =6(.M 


- ] if Ve 
0 if Mev. 


When s<¢, one can consider the simple case where the additive function of the set » can 


where 


be represented by a Lebesgue integral 


w(M.s; v,t f( Ms: Q.t)dQ 3 


where dQ is an element of the »-dimensional space and is a v-tuple integral 


When this assumption is made, one finds, in applving total and compound probability 
theorems that the Brownian displacement function /, which govern the v-dimensional Brownian 
motion,’ satisfies the Chapman-Kolmogoroff functional equation 


{(M,s; P,t) f( Mis; Ow f(Q.u; P.t)dQ ! 
Jt 
for su<t. V denotes a v-dimensional Euclidean space and / is Lebesgue-measurable in V 
In the problem of chain probabilities, {(.M,s; P,t)dP is the elementary probability that the 
state of a physical system, fixed by » parameters, and represented by V/ at the instant s, be- 


comes at the later instant t one of the states of a certain elementary set dP of states. However, 


in this interpretation the transition probability or density probability function /(.\/,s; Pt 
must satisfy the probabilistic conditions 
/ Vs: P.t ) 5 
f(M,s; P,t)dP=1 6 
e/ It 
If the physical system is moving continuously it has to travel the finite distance \/P during 


a finite time t different from Ss. Then if the states V/ and P are two different states we have 


lim {(.M,s; P,t)=0, WAP i 


If the state P coincides with the state M, the transition probability / is no more continuous 


and it tends to infinity as t->s 
lim {(M,s; P,t Ss 


PoM,t-s 


2 By a v-dimensional Brownian motion we uhderstand a 2 mer , ne kofl 1 


Brownian motior 
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Such properties of the transitio probability function when 7 s, resemble those of e Dirac 
function which 1s a kind ol impropel runetion It ollows that the conditior S must be 
excluded from chain probabilities Therefore the ass imption ol the existence Ol al elemen 
rv probability /(\Z,s; PjdP cannot be accepted completely and an exception must be 
made for the case s=t, 14=—P. If fis interpreted as a generalized function, in the sense of the 
distribution theory, then the cas tf and \/=P is covered by the statement that / has 
unction behavior in this case 
Using the additive theory, both N. Wiener [23] and P. Lévy [1S] obtained very important 
results concerning the transition probability hunetion NI. Kae iZzk © Kakutani [14, 15 and 
others in a series ol papers obtain ed properties OL thre Brow! mn motion mn a dimensional 
space and in connection with the Riemann surfaces {n interesting probabilistic discussion 
Ol tle veneral diffusion process Ih one dimension Is also riven by \W Feller 5 Recently 


\L. Kae [13], relving on the properties of Wiener measure in the space of continuous functions, 
showed the disadvantage of Wiener measure theory applied to Brownian motior 

Orn the purely yhvsieal side thie two papers by { hlenbeek and Ornstein 21 ind | hlenbec! 
ind M. ¢ Wang [22] should be mentioned 


3. Method of Solution 


To obt in a veryv venel il solution ol the tunctional eq } = \I Fréchet o6i tried a method 


ol separation ol variables 1 the torm 


nd found that the ecessarv and sufficient condition that the expression 9) should be a solu 
ion ot eq 1) is that or the values of wz such that s u ft the functions A (/ / ind B () u 

rim biorthonormal system ove | Lf 

. 1 ou 

| Ju B (2 wd where 10) 

‘ ‘ ‘ 

) otherwise. 
Nevertheless Fréchet real roc that the expresslol be) does not satisiv the eonditions o 
| 
( probabilities d ro iced ste id the infinite se es 


Ft ge >> A,(M,s) B,(P,t 


mctiol 


f 
a 


nad proved that the necessary ind sufficient conditior | 
\Zs; Pt), satisfving the Chapman-Kolmogoroff eq (4), could be written, in at least one way 
ol f is that the functions | V/s 


the form (11) where the convergence ts uniform 


B Pa orm a comple te orthonorm il set Ove! | 


The Series sol tio}! 1] is s ipposed to bye unilormls convergent lol ( [ fixed Sie 
{ f and VP varving arbitrarily over | 
4. One-Dimensional Case 
( onsidering the cas ol oo limensional Space | ind with thre mual notations 
d or M, P, and &, the Chapman-Kolmogoroff eq (4) can be written in the form 
mi if / / 4 / 4 2 
rechet s method le | is to tuke tor the solution of eq 12 the infinite series 








where the functions A,, B, form a biorthonormal set over V 

The case of the finite interval V(0,27) has been treated by M. Fréchet [6]. The solution 
of eq (12) on (—@,) in the form of a series of products of two Hermite polynomials has been 
given by the author [8, 10], and its analytical and probabilistic properties are discussed in 
detail. 

The most general solution of eq (12) in the interval (0,©) is also given by the author 
[8, 9] in the form 


) 


r+ - _ : : : 
K(2,8; YO = DWH? (TV (y)0"(s,t) exp (- 2 \(ry)2 9), LS (2) LS (yen(s,t), a he Ls 


n=(0 = 


where the infinite sequence of L*-functions {y = (z } 
vy, (2) A exp| 5 fe ‘ed f 


form a complete orthonormal set over (0,©), and L@ (x) is the nth associated Laguerre poly- 
nomial defined by 


In addition, 


1 
A(8,t) :. lor s<f 


a(s) is a positive increasing continuous function #0, chosen arbitrarily, and 


where I denotes the gamma function. 

It has been shown [8, 9] that, for s,f fixed such that s<¢ and z,y varying arbitrarily over 
0,o ), the series solution (13) is absolutely convergent in (0,©) and uniformly convergent in 
(e,~ ), e being an arbitrarily small but fixed positive number 

Using Hardy’s formula 


SO , ™ = ‘ , a Hl 7 y 2 fa J 
3 Ae Li (2) Le (y)O"=(1—8)-*(Oxy)~2 exp o || 1-0 f 


where J, denotes the modified Bessel function of the first kind, as the generating function 


for the series solution (13), we get 


, a | a] } (Ary 
S (ays; y,t a] D afl exp - 21 . / ] ] ! =e } 14 


The modified Bessel function J,(2) is defined by 
I.(z)=17*JFa(% : & 


where z is positive and a is an integer. 

For @ a positive integer or zero, a condition which is not incompatible with the assumption 
a 1, J.(z) is positive. Thus, according to (14), the transition probability / is always positive 
for any values of z and y in the interval (0, Therefore, the solution (13) satisfies the 
condition (5) provided @ takes positive integral or zero values only 


It has been shown [8,9] that, as a consequence of Kogbetliantz’s inequalities [16], the solution 
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for any finite and positive values Ol 2 


}) satisfies the special probabilistic condition 


lim f(x,s; y,t 0, 


and wv, and also the condition 


lim f(z,s; y,t 


5. Limiting Behavior of the Probability 


Let us investigate the asvmptlotic behavior of the probability w 
when the time f tends to infinity To do so the existence ol 


f defined by 


: : Pals , 
fir s: yt exp| , /] ry 4 a SSA & Ty 
) i ad f | 


bility the uniformly convergent series 


should he sought Having inh mind that in the inctor @ 4 ais) a 
positive increasing continuous tunction o tf two cases Mav occul 
] lim a(t 
hen 
lin fist () 
herefore. for x. y. and s fixed and ¢ tending to infinitv we have 
a 
| ; | a ] / CXp J é 
w! I depe as on na y It follows from that the lim 
i he inte 1 (0 IS OL ty 
lim w(2,s; Vt)=22> 1 (S41)1 exp 
? 
pro, ded y 2 This limit cle pel ds on thre initial State whic 
—T Tl is termu olow is due to NI Kré< het, al d it me 
limit 1 the ordi irv Sense when f tends to infinity 
By assumption, a es positive integers or zero values d ther 
Is fulfilled 
2 Che onl Ite ive to the ¢ se | Ss tf 
ait 
\ ‘ 2 fy ad posil Ve | ens 
j ad ] 


follows that 


is 


ld defined by 


s called 


ns that thie probability 


onos 


tends to 


f the limit of the transition proba 


lator \ 











hence 


, , f : r+y , : 
lim w(zr,s; V,t) (ry)? exp | - -—— ] >> Az fay(s) "Li (x) L& (ydy 
Jv 


t>« 
> K, [ai(s) "Wi (x) 
where 
K.= [ve wndy 
Therefore the limit of the probability Ww depends completely on the initial state zat the instant s 
This second case is also a nonoscillatory case. 
The above results can be stated as follows 
Tueorem. If the transition probability f(2,s;y,t) satisfying the Chapman-K Imogoroff eq (12), 
is in the form (13), then the limit of the probability w, when t tends to infinity. exists and de pe nds 
on the initial state . Lt _ the probability LS nonoscillatory, provide d a tak S ro oj posit ( integral 
values. 
6. Stochastic Partial Differential Equations 


A. Einstein [3] showed that the probability density f satisfies a partial differential equation 
of diffusion type whose coefficient is a certain diffusion constant related to the viscosity and 
temperature of the medium in which Brownian motion takes place, and to the radius of the 
Brownian particle (assumed spherical). 

The existence of the probability density f has been established by many authors from 
different viewpoints, especially by A. N. Kolmogoroff [17] and W. Feller [4]. Kolmogoroff 
proved that any function which satisfies the Chapman-Kolmogoroff equation satisfies also a 
second-order partial differential equation of parabolic type, and he obtained a class of part ieulas 
solutions and showed that these, subject to certain special conditions concerning the order of 


the magnitude of moments, satisfy, in variables z and s, the equation 


) o*. f 


of . Of 
or *~ 


+ 0 LS 
Ou Os 


B?(x,s A(x,s 


where 


, ] i 
B?(x,t)=lim (y—ax)*f(ax.t.t t.ydy, 
ui, vo 2At . Y~ ATI At, yay 


A(r,t)=lim (y—zx)/(z,t,t+At,y dy, 
M0 . 

under the assumption that the ratio of the moment of order 3 to that of order 2 tends to zero 
at the limit At=0. 

Kolmogoroff’s partial differential equation was derived by C. M. Techen [20] independently 
of Kolmogoroff’s method. 

It can be shown that the solution (13) is different from those discussed by Kolmogoroff 

To see this, note that if the solution (13) belongs to those of Kolmogoroff, then there would 
exist two functions A and B such that the solution (13) satisfies a partial differential equation 
of type (15). Calculating the derivatives of f, given by (13), we get 


~ r+y @"~* O80 
. a —_ ‘ yy ] a ] @) /»,) 
(ry)2 exp / 4 (J 4 
iS Y | 2 n>0 : A, Of ) y 
O/ ~ r+y 9” a df 
. , a A: ‘ yy _ a) | a ’ a , 
= exp[ 2 lz , ( 279) BO) +L. (2 |e ‘ 
O° / « r+y 6” " a r a’—2a atl 
. . a _ ‘ yy a . a ) 1 ae a " 
dr? (2p)? exp[ 5} a4 L, (x) +(“ r)L J t-( ear. 5 )L r) |L& (y), 


where prime denotes the derivative with respect to z. 
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Kor rand y varying in (0 fixed and siving ina bounded interval such that | 
obtained by differs ntinting 13 ter! by tern are untformly convergent 


ibove derivatives lI LS one vt 


where { r.s.t) denotes 





Since the relation (16) must hold for anv value of y in (0 the « xpression in the square 


brackets must be zero for any positive integer value of n: that is. after dividing bv L‘® (a 
| 


iO y é 
| OO 
\s Is independent ot? and if ~ clear that thre relation iv does not hold inless the 
/ OX 
two coefficients A and J? (4 are both zero Therefore 6(s,f a(s)/a(t) will be inde 
pe rele nt ots which ims contrary to thre assumption that @ is a positive increasing continuous 
function of s This proves the Statement 


It can also be shown by computation that the relation (17) ts possible if the two Tunctions 


Ala and B(r.s) are both zero In fact for 7 0. we have 


On the other hand setting 7 l and 7 2 in relation 17) and eliminating between the 


two relations obtained, one gets 


W cre 

2 r / a Ja r | 21 2 r J J a 2 

f 

i t | Ly 2 ! 2 r 2 i a 2 a 

ibihe 
2 u x 2 2 
D | 
2 t Z / r 2 r r | 2 

The two ratios IS) and 19) are different, hence the relation 17) does not hold unless the 


two functions A and B are both zero 


Now we propose to obtain a derivation of the corresponding parabolic partial differential 


equatiol or the so-called backward equation ol diffusion t| eor\ Kor this purpose we can 


use the generating function (14) and get 


and / satisfies 
>) 


he derivative with respect to 


pl ime denotes 
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Taking the derivatives of (20) with respect to x and computing the different terms of (21 


one gets 


, = (1+6)(r4+y)1f ov it+¢ . 
21{=2r0? (1—@) exp ETy) |i f |, 
2(1—8@) Oz * 2(1—6)° | 
ey eS . 1+-@)(r+y) 24 {/1+8A 2 
=) be 47° (] Q) exp { J Y oO * O ( | l ) f | 0 \ | a 
2(1—@) L Ore Oa\1—@ 23 4(1—@)? ° 4r(1—0 
2 2 = (1+@)(r-+-2 a? (1—6)?+-492r 
(a?+-2?)J,=62 (1 6) exp]. : ” | ; . zy 
2(1 7) ] a] 


. : . ; = * 1+é@)(r+y 
Substituting in (21), and dividing by @2 (1—é@ exp| T “~ | we have 
. { ‘ad 


a 


a OF 4 [r+9) , Of, 21 +6)? +20(1—6)—a*(1—8)*—40ry 
z - 4: 4 7 f—() 29 
ont * 1—@ Ou 1—6)* ne 
The relations 
O/, s o/ 2) 
' and 
Ou Ou Od O# 
, of O/ 
lead us to express f in terms of and as follows 
or O@ 
2(1—86 O/ O/ 
/ ; ( A y ) 
: 1+A)a a(l #)*+YA(1 iW DACs y O# Ou 


Substituting the above expression of / in (22) and setting 


B(x, s)=22? (1—8)|(1+-6) r—a(1—86)? + 20(1—8 20(2r-+-y)| | 27 (1+-0 
27 fa rc a 1Ary 
and 
y(r,8 ri 2jzr(1+é)+(1—0 1+6@) r—a(1—@)*+20 6)—20(r-+y)|\2*(1+86 
2r(1—é a’ (l1—t Ary | 
one gets 
oO? f Of O/ 
6b(2z.8 Vi F,8 +f} () 
Or Ou O# 
On the other hand 
Of Ofd0 a’(s) Of 
Os O# OS att OO 
so that 
OT ais) att oO} a Oo} 
“ 
06 a(t)a’(s) Os a’(s)O 
where prime denotes the derivative with respect to s. Hence equation (23) becomes 
O*f Of. a(s) Of 
8 4 + ). 24 


Ou” dhe Or ' a’(s) OS 

This derivation presupposes certain differentiabilitvy conditions on f and a(s), namely 
f must have continuous partial derivatives of order <2 with respect to z and s, and the right- 
handed derivative a’(s) must exist Therefore the following theorem can be stated 
THEOREM. If the transition probability function f,. defined by 13). is of the class C?, and the 
right-handed derivative a’(s) of the iNCrEaSing continuous function a(s) exists, and is everywhere 
positire and differe nt from zero, then the transition probability { satisfies. in the variables x ane 
Pi) 


; 
? 


s, the parabolic equation ( 
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lt can be observed that the parabolie eq 24) to which the solution (13) satisfies is different 


from equation 15) given bv Kolmogoroff 


oe O/ 
The coefhicient of in (24) suggests setting a(s e, and then putting a(¢ e’, we have 
Os 


@—e 
The special cause d ¢* leads us to a particular solution of ec 12) an the for 
| | i 
firs: yt exp j y)/2\(ry +> " eXp |? / L ta L “"(y 


. [ | 7 ] 7] 7 “9? Ory | 
i 7 CAP /,, 25 
[ (1 iW | 4 | 


und thre eq 24 will be reduced to the normal lor 


O*/ O/ O/ 
Dp ry + () 26 
Ou" Or O% 
If we now set. more cvenerally 
dts ‘ 
where the variable S increases with s, then 
atl f 
and 
A ‘S 


where wGS) is an increasing continuous function of S We then have 


The function f(r.s; y.t) becomes F(r,S; y,7), and the ¢ hapiman Kolmogoroff equation will 
replaced by numbers S and 7 and the function F(x.S: y.7 


still bes tisfie if the instants s, ¢ are 
can be written as follows 
FiaS 7] / kirS { / { / / / S { / 27 


The serie solution 13) now becomes 


yt Pip ‘S JA / I rus y. 1 
J jj : . ' 
exp . | ry pie A, exp [n(S—T)|LIY (LI (y (28) 
When the time t--e, 7 tends either to infinity or (7 being negative as well as S) to zero, 


then @ tends either to zero or « so that the limiting behavior of the probability function @ 


defined by 
Fa S; 4.7 )dy, 


will still be nonoscillatory 
There is nothing arbitrary in F(r.S y,T : defined by 2s which Is a parti ular solution 


of (27 Therefore the most general solution will be in the form 
firs yt fk r.n(s ynit 
where n(s In a(s) is an arbitrary function varying in 
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Some Higher Order Integral Identities With Application 
to Bounding Techniques’ 


J. H. Bramble ' and B. E. Hubbard ° 


July 17, 1961 


Let PR be a simply connected region in Ey with smooth bounding surface S For a suff 
ciently differentiable set of functions we derive a class of quadratic integral identities relatir 
surface integrals of derivatives to integrals over R These identities are a generalizatio ot 
i first order identity given by L. Hérmander (Compt. Rend. Douziéme Congr. des Mathe 
maticiens Secandinaves Tenu a | 1, 1953, pp. 105-115) and L. FE. Payne and H. F. We 
berger (Pacific J. Math 1Y5S) pp. 551-573 \ in example of an application of these 
identities we consider a solution of the boundary value problem Au pu F in R and 
fon s Here A denote the I iplace Operator and USp We obtain point Wist 
prior! bounds for the derivatives of u in # in terms of i quadratic functional of in arbitrary 
functior Hence the Ravleigh-Ritz procedure can be ised to make the error arbitrarily 
ill 
1. Introduction mental solution in Green’s third identity, thus 
introducing more singular behavior of the bounds 
In a recent paper [6]* L. E. Payne and H. F as the boundary is approached In our paper the 


Weinberger give a method for obtaining bounds for bounds for derivatives of all orders have the same 
behavior near the boundary, provided the boundary 
data is sufficiently differentiable 

Extensions of the fundamental identity (2.4) can 
derived for higher orde! operators by the sane 


solutions of second order elliptic boundary value 
problems ly 
see eq 2.4) of an integral identity of Rellich [7] which 
was essential to their method This generalization, be 
which was first obtained and applied to hyperbolic technique used here. One sucl involving the 
operators by L. Hérmander [4], displays the highest biharmonic operator 4’, ts applied in a forthcoming 
derivatives in the form of a second order operator | P&per ol J. H. Bramble and L. E. Payne [2] to obtain 
Lu=(a"u This fact makes the identity useful | explicit) pointwise bounds in the biharmonic Di- 
in treating boundary problems for the corresponding richlet problem 
differential equation a= fF \s has been previously pointed oul in various 
Use was also made of this identity by Hubbard | Places (see e.g. [6]) bounds of the type obtained here 
5] to obtain bounds for membrane eigenvalues by | Ca0 be used in conjunction with a Rayleigh-Ritz 


finite difference methods technique to approximate the unknown solution (or 


In section 2 we give a further generalization of 
Hérmander’s result, which involves higher deriva- 


tives on the boundary and where the highest order -_ A Class of Quadratic Integral Identities 


that paper they give a generalization 


derivative) arbitrarily closely 


terms enter as derivatives of Lu Section 3 gives 

an application of these higher order identities in Lot 

obtaining pointwise bounds for the derivatives of a Liu au... 91) 

function vu. The particular problem treated assumes 

a knowledge of Au pu ih & Pegion Rand u on the be defined on a region FR in #, with boundary C 

boundary C’ where 4 is Laplace operator and p >0O The svmbol indicates covariant differentiation 

in R-+¢ Because of the important physical appli- which coincides with ,, (partial differentiation with 

gris ee cee oe ae “a — respect to a in Cartesian coordinates A repeated 
index indicates summation from 1 to A The sym- 

trate a use of the higher order identity More | metric tensor a(x) is assumed to possess piecewise 


creneral equation could be tre ited With only tech continuous derivatives of orde! \/ l In addition, 


nical modifications Other methods for obtaining let the eigenvalues of the matrix a he bounded 
explicit, pointwise bounds for derivatives in such away from zero and infinity in FP Hence there 
proble ms have been rivel bi J B Diaz oO} and by exists a positive constant. a. sucl t] at for all real 


» ' roe ] , ‘ . ' > 
Pavne and Weinbergei 6} In both cases the nonzero (3 F and all z in 2. 


method viven involves differentiation ot the funda - ds 
I N N 
' nN { a ‘¢ a‘’s a (2.2) 
Air Air I Re Air J i=1 
D ( ( AF 4 s ‘ 
- Md : , : ! Under these conditions L Is said to be uniformly 


elliptic 


261 








Let f(x), A(z), 
of the second order with piecewise continuous 
first derivatives in FP. L. E. Payne and H. F 
Weinberger have developed the following first order 
quadratic identity [6] and used it to obtain pointwise 
bounds in certain boundary problems. It is clear 
that 


(f™au ju ) " (f™a"’) nt, U4 2/"au Bie 
( fa"), UU) +2( fa" u, ,U,,,) 
Zh ab, Ub) m— 2h", » La) (2.3) 


An application of the divergence theorem vields 


f f"a' , 2/ la™™ Ut vp, (ls 2 f" uo, Li u ir 
ev C : JR 
+i i(f"a"’),,—2fi,a"™ }u,,u, lv. (2.4) 


JR 
In an analogous manner we can develop a homo- 
geneous quadratic integral identity involving second 


derivatives of vu on the boundary (. Now 


(PPA ae, 2%, 11) 0 


(fmAtsat) ; ty xt 2h" Ava ue, yu 
2(f* Aa ts, tt: im) "Au, ,,,( Lu 
(f"AVa*') 2(f'A'a") UU 
QL" A (atsu,,), a (2.5 


After applying the divergence theorem we have 


| Af { f*a 2fal™ | u, ,.U, , vs 
+ (f"Ava*) ,—2(f'Ava™) } Ub ude. (2.6) 
JR 


The corresponding identity which involves covariant 
derivatives of order 1/+-1 on the boundary is 


i i rm_ kl o st km . 
p Aj? eee A va fea 2) a UL ki wiht w? wl S 
e c 
‘ pm hij tuju . 
2| : ace (Lu) i, tip... samt 
JR 
mr i,j i ki 
+ (f Ay an a he 
JR 
‘ tl agisj im) km 
2 fA... Ans a )im} Ua... inktip.. jal 
9 m ii) im) kl . 
- J Aj oo LAM (Aietha)ia,. yl; _ 
R 


{ 


) 


~I 


Since the integrands are in each case tensor in- 
variants, we may perform each integration in the 


262 


Aj}(x) be symmetric tensors | 





most advantageous coordinate system. Since bound- 
ary conditions usually ure given im terms of normal 
and tangential derivatives we shall display the deriv 
atives appearing in the boundary integrals in such a 
form Assume (C' to be a surface possessing V+ 1 
continuous derivatives as viven im the parametri 
form 


X G ws : y* S ) l, , \ oN 


The unit normal .Y y satisfies the conditions 


G*..A‘=0, a=l1,. , N—1, 

NOY | 2.49 
and the orientation of Vis taken inward In whit 
follows Greek indices will alwavs range over 1, 

V 1. whereas Latin indices will refer to 1. _N 
We introduce geodesic normal coordinates (y 


y*) in a strip immediately adjacent to the bounding 


surface. The transformation is given by 
y=Giy',... ¥ y X(y', y* 2.10 
The coordinate y* refers to a distance along thi 
inward normal to ( at the point (y a 
The metric tensor has the form 


Yan Gi. 
ee A (2.11 


The Christoffel symbols which involve the norma 
coordinate y* are 


= | Oy 
ae 2 oY 
N p \ 7 
Net ¥ \ — 


The various second order covariant derivatives Ol u 
then have the form 


= Or k in 
Pie Ue. 
oy*O"r ap 
- - Ou Py - 
uf u = . u 
= ie Oy" Oy" aN 


- O71 — 
Wenn ov") P 2.13 
We indicate with a bar when the tensor is expressed 
in geodesic normal coordinates. Our operator takes 


the form 


Lu (a'?u,,) (a U, ) 2.14 


When 4 is the Laplace oOperatol A we | ‘i The following notation wi bye ised throughout 
thiis papel Let fia he nu precew contimuous Tune 
An ul rT Geu 2 15 tion on J?-+-¢ then 
Since ¢ O and « | Also thr Hirst mvariatht PrhaaN 4 
tukes the form ef? +-¢ P efi +-( = 
grad f u c u u 2 % Kor \ 2 ty / re pres i iu Lii¢ closet 
rve Equation » 10) becomes 
We now express thre bound iryv wnitegra ippenarimic 
2.4) in geodesic normal coordinates x GS \ ~~ 9 95 
( % 7 oT > 3 & os nO a where we have set y ‘ urelel i! iiong ¢ and 
P / \ Cistane nulong thie nit imnward normal 
7 = - af. If AGS) os the curvature of C then we have 
a ul 2( feu i / iS 2S ia 
] IK \ \ (} ] 
if L=A tl a nd (2.17) becomes 
2.26 
q*a..i 7 2(F°%..)u.. dS. 7 
r ; Phe ise 0 this coordinate syste is ol course 
218 restricted to a certain strip in the eighborhood of 
the boundary where in terms of it all points are 
he suriace ntegral ! ant Cul 0 tTrented ! ' maiquel\ defined 
Sitnitial manne! to obtain 
3. Pointwise Bounds 
: : di cities ; . As was previously noted, Payne and Weinberget 
! road isc of 24 in order to obtain bounds Ta) 
2( fu a u hs 2.19 SOlutIONS © second order ¢ liptu partia differential 
= equations This entity enabled tiie tO estitnate 
In prartpetitalr it / | ms chosen on € so that thre intevra Ove! i closed Stirtuece Of Utve Square ol 
the hol itdermvative of a solution of a second ordet 
fs 4A , equation ih terms of integrals of the squares Of the 
» 1) funetion and its tangential derivative over. thre 
fa A () a Ssturiaes In orde! to bound tiie derivatives at 
pol t within the region the oted that one could 
then tor / A 2? 149 tukes the su ple ol ameres Liste Live “Green's Ider lil bye mie « ireful 
to define a parametrix im such a Wav that this is 
° permiissibi The resulting esti ites have the dis- 
u ul —_ iS , idvantage that the coefficients become infinite more 
. rapidly ns thre boundar\ is ipproa hed than do 
those i the estimates for the value of the solution 
In view of Z.15 we have therebv isolated an u itself This procedure has also beet! ised by Dinz I3 
teresting combination of mixed normal and tangen in estimating derivatives at a pomt for solutions of 
tial second derivatives on C' in terms of second order Laplace’s equation 
tangential derivatives on ¢ and certau Integrats Qne could proceed mann alternate anne That 
ove! RP is, Instead oO differentiating tive Cireens ldentits 
In fact, a specific mixed derivative, say Wye, ca simply write the identity using for the funetion to 
be isolated in this manne Let be evaluated at the port the derivative of the 
solution This leads to the estimation of surface 
fj | ji {) integrais OF squares Ol second derivatives of the 


— — function ove! the surface in terns OF Ssuriace mite rrals 


{) — »>o 
A I —— ol squares ol tangential derivatives ol thre hunction 
over the surtuce The identits 24 vives a means 


‘her | 4 + | ' Tn ~ | 
en rol i \ (2.1 bec : of obtaining the necessar\ Inequalities 


° As an exan pie oO; an application of these identities 

= 2 fe. 3120 292 we consider the problem of obtaining a_ priori 

4 pomtwise bounds for the derivatives of a function 

i for which Lu=Au —pu is know! R+C' and w is 

As we shall see in the next section, these choices o known on Here R is a simply connected finite 
j A will enable us to obtain new pomntwise bounds plane region bounded by the smoot! closed eurve 
for derivatives in certain boundary problems ( A s the Laplace operator md plryy >) and 
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bounded together with its first derivatives in R+C. |) Now on (C, Au—pu—uwy+ viss—pu, so that " 
(The function p is taken to be positive in R+C for 
the sake of convenience. By a slight modification, 


similar results could be obtained for p>0.) We b MMi AS = Au — pu " 
assume that Lu and wu are sufficiently smooth in | ~ : I 
R and on C respectively. 
As a starting point we use the “Green's Identity” WoalS +2 uredS S 
I 
(P) PAvdA-4 | (v2 —r oe )ds, 31) | Noting that 
JR - ( on on 
: a Oru - Ou OU - Ou 
where v 18 any sufficiently smooth function in R+ C, u ds A aN ds? A a) J 
P isa point in the interior of 2,T 5~ lnrpe(? re= dis- 
-7 we have that 
tance from the point 7? to another point @), and ' { Oru ou) 
on u u“ 1a<2 | || pu > Au pu ) 4 ( > ) j 
is the outward normal derivative. Now let | °° si \ 
ou ou\? , a ’ , th 
= . Tv " " " T ‘ ama / iS > | , 
v U, => Then we have thy Vf ( “s ) ! 2g l 0) k 
ae “ or On, It follows from Green’s Identity and the arithmeti 
will) J ram “f Mot on on )ds, (3.2) reometric inean inequality 
which may be written Vent b ear b ; 0) 1] : 
x, (P)=u,,(P) P(Au—pu),dA (pu), thiaat 
JR Jk ) 
, | > ow : 
or . du. ypurtA+D(u.u pu IS 
+ | ( “ yds. ? 4 Ji 
Je on On . 
1(° (Au— pu ry ; 
We shall consider x,(/’), sinee x,(2?)—u,,(7’) is 2 | ' 
assumed known. Making use of Schwarz’s inequality 
for vectors it follows that Combining (3.10 3.12) with (3.4) we obtain 
me . ; 
xX (P)x,(P)< kK, ( = peda \-D(uu x, (p)x,(p): h ; Au—pu HA 
“JR <-J 1 P 
: . “rf . 
Ou., Ou ; | Oru Ou ; Ou 
Deu,,U. dS 4 dS 24 {+ ‘? | pu + diN-+ iS 
t q r =e = ) ¥. || ~ ) +{ ds ; ( Os ) 
eee ae lf ll 4/0 ‘ he 
where D(u,u) is Dirichlet integral, and | wdS4 ( - ) dS+2 u IS) 7 
2 2 %\o 
K s| I? grad yp "1d A+ pd A | 
JR JR 


All terms on the right hand side of 3.13 except 
+4 (2) as4 p rads. 2 5) the last two are in terms of data In order to bound 


Cc on the integral of the square Ol the normal derivative 
around the boundary we write the identity 


Now we have 


i Ou Ou yp, OUOU] , Ie 
f Ou,, OU, . I (( 5, )—{ mn) i Oo) fas 


dS < u,,u,, dS 3.6 “ Wi 
on On J « . ‘ vect 
oc 
a. / {7 \u,,u, dal 2| fu, AudA 
. . . CTOs 
and in terms of the normal coordinate system intro- vt vt 
duced in section 2 (with the boundary coordinate 5.14 : 
taken to be are length, Le., y’=S, Y=N) 
° - which has its left hand side essentially in the for 
. : 7 9 — ‘ ow e. are . 
uM g@S ivy + 2ut, 28 (3.7 of (2.18 Here a i ny and are the unit nor 
Je Je mal and unit tangent vectors respectively and /* ts 
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e bounded together 


an arbitrary continuous, piecewise continuously dif We assume also that fo and AY as 

ferentinble vector field in R4-¢ We choose /* such | with their derivatives in B+ ¢ If the boundary 
that /*) on ¢ For example, if R is star-sh iped | is sufficiently smooth one can define these tensors 
vill respect to the Origin then we mav take / J ina thee SUA vstem) as follows Let 


Krom ,.14 wi 


where ¢ ims any 

the coeflicient n 

K ire positive 
K,+-A2=-5 


/ 


/ 


nihiedial 
pb (S")as 


Kin il \ Wwe Tie 


Wi i ike lise of 


vector field f' ar 


obtain thr Inequality 


| 
for N<s 
/ 


4( /“s Ou 
fr, 1(2") LS j | i) {or \ i 


. ; 
AK) D (uu J nud A - 
A JR where A is a constant greater than the maximum 
. curvature Ay of ¢ 
' ; pJ'(Au— puyida 15 Using (3.11) and Schwarz’s inequality we obtain 
A alla irom SIS 


bound ior the leat 


utrin ("0 / / and AY and | 7] LS | LS + Au— pu Au— pu JA 


est eigenvalue of 


Let « mia if iP) and choose ° . vs 
For example take A . 5 - t | vu pu LA | « t LpdlA (3.19) 
‘hen we have. using (3.12) and (3.11 where 
( 2} "Ag "A 2( f'A 20 
, | E | The last term on the right hand side of 3.19) 1s 


bounded is follows: 
( Hat AACKBY u,,,u,,aA 3.21) 


where Bisa cor stant After a moderate « alculation, 
Making use oO the definitions of j and A in the normal 
coordinate systen if is possible to obtain the bound 


B=1+K—K max K,2 2(K’ Kk. K |r) 


p au where A’ is the derivative of the curvature with 
respect to are length Now usInY the divergence 
&, theorem we have 
a) 
. ) 
| //s l ‘ Pa u, (x n )u,as 
| JR Je Ou ] 
//) 


where hie } are Cartesian components of the 


od to obt mn bo ind lor ne* availa ® 
‘ exterior unit normal We may rewrite the boundary 


2.23) and (2.6) assuming that the integral in terms of normal and tangential derivatives 
d the te sor field A have been ind obtan 


chosen as indicated We tl Cl obtain the identity 


~ Ou O7u 
uu. aA 2q : dS 
R 


Z ("Ady Au iA ; Je On Os? 
gi—2(f'A Ui ,u, aA 18 | p K ( ( > ) 1 ( + ) ads | } (Au)?*dA 5.23 
J . . JR 








Using (3.11) we have and hence, using (3.17) we have 


Jk 


. ’ wae . % - 3? . . ° | 
|, uy MM ASP K =) dS + p ( = ) dS pu (pu AA<¢ | dS ” 


‘ F ou? ... ; ‘ . 
! +P (1+K) (Si) ast2 J (am puyda re po") as es) (an—puy'ds (3.29 7 


, ‘ . where 
r2 (puyrdA. (3.24) Va 
A 1) ) 


Using (3.11) and (3.12) we have the inequality ha 


(mustd<pu 4 wdS - p Mu 
JR -'¢ und 


+q (S ) dS4 ( = _— iA}. (3.25 2p TE IE pw 


‘ 


° . = 
, , = , Finally we have from (3.9 3.11), and (3.17 pr 
Combining (3.17), (3.24) and (3.25) we have } 
‘ . , 4 ' . "(Ou i 
4 Z > 2 weeds <2K 2a q wdS IK a dS 
| u, ju, AA< b Putas +b | ( on )ds f . Ji " | ( Os ) 
K Jt Jc \Os 
* ty \? . 
(Ou? ,. : 1.2 | dS +2Kh4,a Au—puyrdaA 3.30 
h (> : )ds ' b, | Au pu A -) ( OS ) J R f 
Jf _* JR S 
ssieore Combining (3.19 3.21), (3.26), 2) mine 3.30 
— we have x 
h, { q Dnt 2 ns 4 | ; 2 doe 0 » . . 
, 7 aes Kong = | 7: 7 iS IS / =| AS ~ ) S 
+ Ky, bg =1 and by=24+-[(1-+ Ada + 2palay4-2pu/p mm | ( Ox ( Os 
3.20 
‘ Au puyrdA (A\M pu Au pu | 
Now in (3.19) we are left to consider the terms Ji Jt Sin 
% F j Whi re Sel 
WedS and pu), (pu) dA, sio 
d 2A yl Bh ( f oh yl Bh ( 
In order to bound the latter we can write ; ; 
d y 4 Bh ls 2A wil Bb, c, and d 
| (pw) (pu) dA Now (3.17) and (3.31) mav be inserted into (3.13 


to vield the desired bound wh 


(p.,p.,u? +2 pu. p, ut pu.,u. dA F a ) ord 
|, / pP / pP / xX P F K {6 wdS 8 | ( Oo / ) IS 
e Jc\O8s tity 
: [ (72: P, Lp )( tu? en VA 497 * />tu . 
JR pP 2 b ( ds? ) dS +B Au pu AA 
S eo R 
op p 
where (3.11) has been used with a =P i Pst, . > 
3 8, | (Au—pu),.(Au— pu), WA} (3.82 
J R 
Now (3.27) and (3.12) vield No 
where the 8,'s are explicitly determined constants 
lt should be winted out that lor tl e sul ee oOo 
(pu)..(pu) dA | ! \ ) 
|, - 2 sunplicits no effort was made to obtain the “best J i 
] 2p, P, ; : _— “/ou\? such bound Better bounds could be obtained. for 
5 ( > P ){ f wd + t ( on ds example, by leaving certain known quantities under 


the sign of integration (compare equations (3.16 


| { Sa pu)’ j i} (2 98) and (3.17)) rather than replacing them by a mani Thi 
de 5 teri 
p 


JR mum Value 
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4. Higher Derivatives in N Dimensions The boundar\ rite erals Occurring i the ibove ili- 
equalities are the same as those to be dealt with 


TI c proc Css des ribed im section 4 cn he vreneral below 
ized easily to higher derivatives. In bounding the We shall show how the desired bounds may be 


Mth derivative we must require the existence and obtained for N<3 since in higher dimensions only 
] . 

integral ilityv of Vth derivatives of Au pu, pu and technical details entel F 

\/ 1 tangential derivatives of uw on ¢ We also Applying Schwarz s inequality, as was done in 

issume that bounds hav bee! obtained for the 3.4 we see that 


various lower order differential invariants 


The mode of reneralization will become upparent x P \ P K { pu), pu), 
rom the case \J—2 where we assume that bounds L JR 


nave been obtained lor 


‘ Ou : — ray) 
purd Dau u tal dS Uo U AS ~ ( ~~ Or J 
. . . QO! > 
1] where Afr is given by (3.5 
We see that 
is Was done eXpil tly ! iwo din enspons thre - 
previous section Du , Mia (4 P, ) pu 
| ; pr pr j ‘ ] 
1! denotes the surface area of the unit sphere Ji Pp ud R 


V dimensions then 


but by assumption each integral on the right side of 


San indamental solution of Au—0O We ca write 15) has a known bound Sinilarly we assume that 
' . a bound is known for Q@u,,,u,,,dS 
) , ) j a 
Xis(l l (au ' | | on It now remains for us to obtain bounds for the last 


term ol 1.4 We see that 


OI Oi 


Siice l / does not enist for \ -towe cannot apply % = = . om = 
J R | greg? .+wgs ndS 2q GPU ann) prwd S 
Schwarz’s inequality, as was done in (3.4), for dimen- ‘ : 
sions higher than three ' 
For w\ } we appl Schwarz’s inequality as follows = aS 16 
( I ul ): / I d / und From 2 15) we see that 
ii ew = (AU— PU) e+ (pu 7 7 
where the singularity of T is chosen as the origu It 
| The first two expressions on the right are known and 
order Lo hound j und we use (;reen s Iden- 
titv”’ . of ou h - 
u - u 
OY oY" OY ad 
4 Awe I k’) 
. | ha 3 : te , (4 s 
a Dp’ 
Boundary Integrals ’ 
| where 
Now us « 4 | we can ensil\y obta n Ou 71 l j 7 . 9) 
OY yh , 
} 
Ta 7 ) } _ . , 
; XN ri , - Pp? Ihe first term on the rigl { side of 1S) is known We 
can use the inequalities 
Boundary Integrals. 
hed, U4 < 7, gu, a 
‘his Sibiiic basi technique tite be used to handle all 
terms in (4.3) for N>4 bP Mia < Tog gh aM) j2, $.10) 
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where 7;, 7. are upper bounds for the largest ¢ 
value in each case, to complete the bound for 


fp JU annUayndS 
Jc 


by a function of the quantities (4.1). 
Furthermore, from (2.15) we see that 


Urn (AU pu) nT (pu v— JU aan 
and hence we can bound the last term of (4 


terms of known quantities and a bound for 


f g® gi, avy sand S. 
c 


From (2.7) for 14=2 and the choice 


{* ] J x 0) 

Ab g”, Ai : AS N 0) 

Be’ —1. Bv=—o i7J#N, 
we have 
fp [98g U ayn pan— JU ann Uavw dS 
ef 


— f"A i 3 | (Au pu) + pu 7 rm 
J R 


+ (f" AB 5"), .—2( f° AB), } Wy ipl 
R 


Hence 
p g® Q?Tayv UpsvdS < p g* Ti any UisvvdS 
+ ( AB (Au— pu), a (Au— pu),, 
LA 
+ (pu); a (pu) ;jdo+T(u;,), 
where 
T(uii4) = / { (fA BY 5), —2(fAvB),, 
JR 


+P PAT B™ } Wy ap Uj rsd. 


1gen- 


(4.11) 


6) in 
$1.12 
dV 


(4.15 


Sut 


where 


Is 


again 


eigenvalue of the 


~ 


ul 


coefficient 


di 
“ ) 
yils Au f do f 
u,.ds—r Au tia 
o/ R 
+ € Gg i rT ds } lo 
1 upper bound for the largest 


matrix 


inh £.15 Wi 


can substitute (4.1) into (4.4) for « l and aches 
the desired bound for 
Gg i yll / 
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A Priori Bounds in the First Boundary Value Problem 
in Elasticity 
J. H. Bramble ' and L. E. Payne 


Au 9 Log 
‘ Lear f elast \ Phe i 
i i I 4 ] d e the | i 
1. Introduction 

I recent papel 2 the authors gave a method for obtaining bounds in the Dirichlet 
problem for second-order elliptic partial differential equations. Pointwise bounds were ob- 
tained which involved, 1 addition to the approximating Tunction only undifferentiated data 
OL thre problem The Inequalities were Ot nu priori type henes only regularity conditions were 
imposed on thre ipproximating tunetLor 

lt this papel we obtain bounds for the stram energy and pointwise bounds for Like solution 
ol the first bout dary value problem mi tine equatllor s of elasticity Other methods for obtain- 
l bounds which are based o1 known variational principles I ive been given 1n the | erature 
see, lor instance Prager and Synge [9 Diaz and Greenberg [3 svnge 10 ind others 


\ method based on a priori inequalities was proposed by Payne and Weimberger [8 
Their method involved the introduction of a somewhat complicated unit-tangent vector, and 
pointwise bounds involved differentiated data. We derive a method which vields pointwise 
bounds which involve the undiffer ntinted data of the problem 


1 tiie problem l det consid 


ition the displacement vector with components u,; satisfies 


the following system of second-order equations in some N-dimensional bounded domain J) 
| 


/ f / vil / |] 
Where r involves the elasti eonstunts A ind u a \ uu uu | 2a r denoting Poisso! 3 
ratio ind F Is proportional to the body force In (1.1 denotes O Or,Ou and a repeated in 
dex ndicates summation from 1 to .N. where N denotes the1 imbet of dimensions By special 
izing to two or three dimensions we ecun treat the two physi illy interesting cases In addition 
to thre differential equation ef thre displacement components satis the conditior 
if / | ) 


on the boundary (C' of DD Regions whic extend to infinitv may also be treated | suitable 
conditions are imposed at infinity 


Tl e strain enere’sy I wou lor problem 7 . 2 Ls defined us 


ry 
Qu He uae E l u u u a (u qa )) 


ur problem then is that of obtaining uppel and lowe! bounds fol hk UU and ior the value 
of u,(P) at a point P in DP. As we shall see, these bounds follow immediately once we have 
established eertain a priori bounds Let us Suppose that we are able to establish the following 


inequality for sufficiently smooth functions w 


E(w.w)<K nw wds+h | n. grad.w,-gradwdS+kK LiAw)L 


. . . 


widV 1.4 
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where », and 7 are given positive functions on (’, and A,, Ay, and K, are explicitly determined 
Here grad, denotes the tangential projection of the gradient. Now if we set 


the right hand side of (1.4) will involve only the data of (1.1) and (1.2) and the approximating 
function ¢;. We would then have an a priori bound for E(u—¢,u—) which could be made 
small if the functions ¢,; could be chosen to approximate the data closely in the mean square 
sense. But by the triangle inequality 


vE(¢,6)—vE(u—¢o,u—@) <y E(u.u) << F:'(¢,6) + E(u—o.u—o 1.6 


Hence knowledge of the a priori inequality (1.4) would give upper and lower bounds for E(u,u 
The inequality (1.4) with w,—u,— ¢, is suitable for application of the Ravleigh-Ritz technique 
for improving the bounds. 

The pointwise bounds will also follow from appropriate a priori inequalities. We concern 
ourselves then in this paper with the establishment of the necessary a priori inequalities 


2. Auxiliary Inequalities 


In this section we employ some generalized Green’s identities to derive certain inequalities 
which will be used in sec. 3. We consider first a function h harmonic in 7). As indicated in 
[8, eq 2.4] the function / satisfies the following identity: 


D [find h of iT] h A dS {* 4 2} h h AV 9 


. vi 


where the /* are components of a piecewise continuous vector field and 6, denotes the Kro 
necker delta. The n,; denote the components of the unit normal directed outward from )) on 
C’. In [8, eq 2.4] we have taken a” =6, and used the fact that 4 is harmonic. By breaking up 
the derivatives of A in the boundary integral into normal and tangential components we obtain 


f f*n, | grad, h 2_( oh | dS ( f*, 6 of jh h dV 2p f' vrad, A-grad, a On aS 2.2 


on J JD n 


oh : ;, 
where or denotes the normal derivative of h on CC. If we now choose the f* in such a way 
v 


that t=/*n,>0 on C, we obtain the inequality 


4 oh 2 — “4 e . ‘ : ° F oh . a9 
P t(5,,) a8: pe grad, h 21S 24 fi orad, A-orad, a on dAS+aDth.h). 23 

In (2.3) we have used the notation 
D(h,h) hh dV 24 

e/ D 
for the Dirichlet integral, and the inequality 

[2f" f*,.5,\h ah. dV <aD(h,h). (2.5 

JD 


The constant a@ is clearly any upper bound for the largest eigenvalue of the coefficient matrix 
2f';—f*.5y. Then from Green’s first identity, i.e., 


Dh, h) | h * dS, 2.6 
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Oh , . 


. Yh | . 
¢{ = t grad, h?dS | 2f' orad, h- gra ] 
O 


Now iron Schwarz ! 
} . . . 5 
, ' : Oh 
f h- ora ahi?dS - | t ( )dS , 28 
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0 rhe rit ao etri ( equ |  ¢ ( equa t\ 
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geometric mean inequality and the fact that @ is positive, we have then 


' Spa : , : , 
f tiv Yi ro, WsldSs| b+ ; ] Vi Mis poh Ws j\dV +8] LAWL Md 2.14) 
Cc { JD , 


for any positive 8. Now since y, vanishes on (' we have 





Wy, Vv tay, yp dV WiLL pdV 2.15 
/ D * 
or by Schwarz’s inequality 
dA tod. &, 1V LAWL. (wal 
Jl ‘ * » 16 
| vw te ee 
Ji e 
However, for a >0 
Ce eee ee Be 
Jl n> ) 7 
vavlV Ful 





where \, is the first eigenvalue in the fixed membrane problem for J). Lower bounds for 


are easily derivable from the Faber-Krahn inequality [5], [6], or from monotony principles 


Thus from (2.16) ard (2.17) we have 
: 7 eT i aed 
| vad: Vitis bows. Wy s]dV: 
JD \ 
The insertion of (2.18) into (2.14) then vields 
p tly Y, tad, Vv, dS <[b/r, +241 £44 \ LAWL AWdV 219 
where we have made the optimal choice for 8. 


3. A Priori Inequalities 


Let w, have piecewise continuous second derivatives in )). We seek an a priori bound 


for E(ww). To this end we decompose w, as 
u Ww u | 
where 
Lu ) 4 
in D and 
Ww u : 
on C. The function w; then satisfies 
Liu Lu } 
in D and 
u 0) 5 
on ©. Then clearly 
E(w,u | OTST: E(wiu 6 
since # and w’ are orthogonal in /-norm (i.e., E(%.w’)—0 We derive first a bound for 
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ke uf uy Fron the tact that u Oon € we have tron 2.138), with y¥ replaced b 


Just as one establishes the Dirichlet principle for harmonic functions it is possible in 


bounding /{(7,7%) to show that 


New) <RUTH 8 


vhere /7 is anv vector which satisfies the same boundary values as w on ¢ Let J/, be the 


components of thre harmonic vector whicl takes on the Sale boundary, Values as uw Chen 


oO ! DOSILIVeE 7 Now sin the 77, are harmonic we avy ust equality (2 


HT, OH 
s<rly 
) ) 


Which Ca clear ay put mn the tol l | } i so des re 


for the u tevral ol the sum ot the squares ol the displacement 


bound 
sw, over D) This Inequality will be used in se > to derive pointwise bounds fol 


Wi seek then a bound for ) V 


lo this end we introduce the auxiliary vector y wit 








By Schwarz’s inequality for vectors we have 


. 


wwdV= 
/ D 


L,OOL wavy pt lap w (wn dS {|x 
/ D vv C -) ( 


‘\ 2 ‘ 
Fax, X%,sdS > 4 | xx. 


Since x; vanishes on C we may use (2.18) and 


| | mw wav P <q h X 2(1 ¢*f ree | / wow 
a/ D ( 


2.19 


with 


+ Ww { 

J \ 

If we should now set w u Qo where uw, is the solution ol l 
mating function then the right hand side of 5.17 would 


We employ this inequality in the next section for obtaining pointwise bounds for 


4. Pointwise Bounds 


In order to obtain pomntwise bounds we as 
bounds for #, are then derived by making use 


homogeneous equations of elasticity Such 


dimensional case by \quaro !] and Svnge [10 
obtained by Diaz and Payne [4] 

Let S(r) be a sphere of radius r and surface > > 
wise bounds are sought. We choose the origin of 


and require that S(r) be contained in D) 


Mean value 


The 


nueaih dee oOlipose 


of a mean-value 


with cent 


In what follows whenever a Greek index is repeated the si 


be used, L.e., Ue Ue will denote the square of the 


ponent Ug. 


From the mean-value theorem 14, eq o.9 


If we now integrate +.2) with respect to? from the 
still contained in J), we obtain 


N(N+2 (" 
2(N+a)R* “Ws olla 


. 


we(0)= 


An application of Schwarz’s inequality then vields 


= “yi? N N+2) PF V(2—a 2a(2 
Wal sale > 
A } 2 N ce - 


A(N)R | 7 @dV 


vradient 


By our construction we see that ? is actually the distance fron 


274 


VL w)Lw v } 


replaced by Xx 


J 
- 


theore 


-. and @, | 


involve / only throw 


i spl rical coordi ite SsVste 


to obtan 





Sal ippron 


e 
| 


ro ~ 
theoren lol the sol ition o| 

nm) was given in the thre 
inalogous \ a Denso! | res lt Wis 
er nat thre pomtat w pot 

| pol 
nitnation conventiol — ( ice) 

i a singel displace Col 


Now fro La since / 0 it follows that 


7(0 \N)AGR : Pe, +-(n)2\dS 


LUN)ACR pt | S 14 


In order to obtain bounds for O) we troduce the Kelvin fundamental solutio: OW 
whi denotes the components of displacement the Ath direetior corresponding to a point 
oad © orig he 8th direction ( 





1) 1) ) 1M) (| e aerivative o ( nctllol / tlie erpre 0 
\ bye ) decon pose ~ Tha poOMLWIS i nas tor 
cle eS « ‘ O ol db se 0 ‘ vtune nequalities \ ( 
oO oO tt ‘ Oo deriy ves oO |) nd Pa } t.14 
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Upon integrating with respect to r as before from 0 to R we obtain 
uN(N-+4) 


N+i4g= , ‘ > . . 
ii ‘v3, (0) ar > lCN +-2)(1 ar-( 230 r,O3R 
Wy( N+-2a+2), 


S(R +a( N+2)(N+-4) aga, 2a ta— N—2)63,27,|0,.dV 1.1 


where S(R) again is required to lie interior to DD An application of Schwarz’s inequality 
then vields 


t3,(0)/?<Kg,(N)R-“ T,w.dV 1.14 


where Kg,(N) Is an easily computable function of .N 


To obtain a bound for 74,(0) we differentiate (4.7) to obtain 
T 0) T Ww, ) aw } IS T L p” lV $15 


where we have used the notation 


af 025(P,.Y) , O25(P.Q 027(P.) aT 


Ou Ou Ou 


the differentiation being taken with respect to the variables of the point ?. Thus thi 
occurring in (4.15) 1s formed by first taking , P (J the displac ement at ( due to a point load 
at P performing the indicated differentiation with respect to the variables of the pot /. 


and then evaluating at the origin. As before, from Sehwarz’s inequality we obtain 


T4(0) 2 ( { t T a" |) ; T n dS 
The insertion of (2.19) into (4.17) then leads to the result 


r’,. (0)|2: { a, N iz w) wa | [ L.(w)dV 


where Ma,(N) Is an easily computed constant The desired bound or 7 1) now follows 


t.1S 


WHY 


since 


T34(O) "<2; | r3,(0 T 0) ‘ 1.19 
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